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1. Introduction 

One of the highlights of the previous Durham Workshop on HERA physics in March 
1993 was the presentation [Q, for the first time, of the measurements of the proton 
structure function F2(x, Q^), which showed a dramatic rise with decreasing x. They 
immediately ruled out a Regge-based description, which is so successful for hadronic 
and photoproduction {Q^ — 0) cross sections. The Regge model predicts a much 
slower rise with decreasing x than was observed in deep inelastic measurements of 
F2. The observation generated great activity at the 1993 Workshop, particularly as 
application of the BFKL equation had anticipated a singular x~^ growth with 
decreasing x. To quote from the Introduction to the Proceedings 

"The argument over the interpretation of the rise began immediately and 
continued all week — was it evidence for the singular BFKL behaviour of the 
gluon; was it just evidence of the need for a different input parametrization; 
what did it imply for the Pomeron and diffractive scattering? Although the 
protagonists tried very hard, it will take a lot more data and quite a few 
more Workshops to answer all these exciting questions!" 

Now two and a half years on, at the time of this Workshop (September 1995), the 
measurements by HI and ZEUS have indeed improved remarkably. Our knowledge 
of diffractive scattering has blossomed and it now has its own Working Group. Both 
the precision and the kinematic range of the HERA measurements of F2 have greatly 
increased. We now have measurements for x < 10~^ and ^ 1.5 GeV^. In fact one 
of the tasks undertaken at the Workshop was a study of the consistency between all the 
measurements of F2 (see section 2). The present surprise is that the strong rise of F2 
with decreasing x appears to persist down to the lowest observed Q"^ (Q^ ~ 1.5 GeV^) 
and yet the photoproduction measurements show that the slow rise of the conventional 
Regge description works well at = 0. The transition from the "soft" to the "hard" 
regime appears to be very rapid. It is appropriate to review our knowledge of the 
low behaviour of F2] section 3 contains a summary of our discussions. In fact 
sections 2 and 3 contain an overview of all deep inelastic data. The high precision 
measurements of the fixed-target experiments are now being extended, for F2{ep) at 
least, to much lower x by the experiments at HERA. Specially compiled plots are 
presented to demonstrate the complementary nature of these experiments, as well as 
to indicate their relative errors. Also we summarize the information obtained about 
nuclear shadowing. 

Partons satisfying the conventional GLAP(Altarelli-Parisi) Q evolution seem to 
be well able to describe the new HERA and fixed-target deep inelastic data down 
to as low as ~ 1 GeV^. Could this be an indication of the precocious onset 
of perturbative QCD and the absence of higher twists? Such a conclusion would be 
premature. At larger x {x }^ 0.01) there are many different types of high-precision 
fixed-target constraints on the individual parton distributions and, with the possible 
exception of the gluon, the partons are well determined. On the other hand at small x 
(x ^ 10"'^) we have, so far, only one type of structure function measurement, namely 
F2{x,Q^), and there are many different partonic descriptions at small x, particularly 
as we have to supply the (non-perturbative) partonic input at some scale for the 
GLAP leading logQ^ limit, or at some Xq for the BFKL leading log(l/a;) limit, of the 
perturbative QCD evolution. At small x the dominant parton is the gluon and the 
description of the observables is driven by the behaviour of the gluon distribution. 
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According to perturbative QCD we expect the small x behaviour of the (sea) quark 
and gluon distributions to be strongly correlated due to the g qq transition. In 
section 4 we overview the different QCD contributions to F2 at small x. We see 
that F2 is too inclusive a measurement to be easily able to distinguish the dominant 
underlying perturbative QCD component at small x. Thus some of the questions of 
1993 remain, but we now have a much clearer idea why the search for the answers is 
so difficult. 

Unlike the conventional CLAP gluon distribution, the BFKL distribution f{x, k"^) 
is unintegrated over the transverse momentum /ct of the gluon. The BFKL gluon 
has two main characteristics. It has the familiar x~^ growth with decreasing x, but 
accompanied by a diffusion in — typically a Gaussian distribution in log fc-^, about 
some initial transverse momentum which broadens as •\/log(l/a;) as x decreases. The 
inclusive measurement F2 integrates over this characteristic k\. dependence and, worse, 
they sample contributions arising from diffusion to the (non-perturbative) low kr 
region (see section 4). This suggests the value of studying deep inelastic final state 
processes which simultaneously expose both the characteristic x and dependence 
of the BFKL gluon. Section 5 discusses two such measurements, the energy flow in 
the central region and the production of forward jets. Whereas the former one seems 
to be very sensitive to (unknown) hadronization effects and therefore, in its present 
form, cannot be considered to be a reliable test for BFKL dynamics, the second 
approach looks more promising: data of the 1993 HERA run are confronted with both 
the analytical BFKL prediction and fixed-order matrix element calculations, and the 
agreement with the BFKL calculation is quite encouraging. 

We devote section 6 to charm production at HERA, processes which are directly 
driven by the gluon. Both open charm, for instance the sizeable charm component 
F2 of F2, and J/ij) production processes offer valuable information on the small x 
behaviour of the gluon distribution. We discuss the status of fixed-order calculations 
of F2, present numerical predictions and try to answer in some detail to what extent 
the measurement of the charm component of F2 can be used to pin down the gluon 
structure function. Particular attention is also given to J/ip production which may 
be able to discriminate between different gluon distributions. 

Last, but not least, in section 7 we summarize our discussions on the spin structure 
of the proton. Many experiments, prompted by the original surprising result on the 
proton spin structure function by EMC, are now yielding much information. After 
a brief introduction of the general background we give an overview of the present 
experimental situation and the status of the Bjorken and Ellis-Jaffe sum rules. The 
small X behaviour of g\ seems to be of particular interest: evaluation of the first 
moment Pi requires an appreciable extrapolation of gi, while recent QCD predictions 
indicate a rather strong rise of gi in the small x region. Finally section 8 contains 
some brief conclusions. 
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2. Overview of structure function data 

2.1. Introduction 

Here we review the status of the data on structure functions and derived quantities. 
Until a few years ago this knowledge came entirely from fixed-target experiments. 
Now it is being complemented and extended by the results from the HERA ep collider, 
especially in the region of low x where the dynamics of a large number of confined 
partons has to be understood, as well as at very high scales where the perturbative 
assumptions can be further tested. 




Figure 2.1. The kinematical region covered by the HERA and fixed-target 
experiments. 

At this Workshop the HERA experiments presented results from the 1994 data 
taking period 1^. During this period HERA colhded 27.5 GeV positrons on 820 GeV 
protons, as opposed to the 26.7 GeV electrons in 1992 and 1993. The centre-of-mass 
energy of the ep collision is 300 GeV. The new F2 results extend to larger values due 
to a tenfold increase in statistics, compared to the data collected in 1993 (and about a 
factor of hundred increase with respect to the 1992 data, on which the results shown 
at the previous Durham workshop were based). Furthermore both the HI and ZEUS 
experiments have made a special effort to obtain measurements at lower values of 
so as to explore the region towards ^ 0. Values of down to about 1 GeV^ have 
been reached as a result of detector upgrades, operating HERA in a different collision 
mode called shifted interaction vertex mode, and by using events with initial state 
QED radiation from the incoming lepton. At the time of the Workshop all data were 
still preliminary, but by now have become final to a large extent ^. For the QCD 
interpretation of the data, such as the extraction of the gluon density, the 1993 data 
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will be used. 

In the HERA experiments, the scattered electron as well as the hadronic final 
state are measured. Apart from tantalizing questions on the dynamics of the hadronic 
final state itself (see Section 5) this allows us to determine the kinematics both from 
the scattered electron and from the hadronic final state. In practice, for the latter, 
a mixture of the hadron and electron information is used, rather than exclusively 
hadronic final state quantities. 

The kinematic plane covered by HERA and the fixed-target measurements is shown 



in Fig. 2.1. Generally measurements at HERA can reach Q (x) values two order s of 



magnitude larger (smaller) than those reached by fixed-target experiments. Fig. 2.1 



shows that the two regions make contact and thus the continuity and normalization 
of the data can be checked. New upgrades of the HERA detectors will allow the 
exploration of even lower in the future. 

The (unpolarised) fixed-target deep-inelastic scattering programme will come to 
an end in 1996 except for the continuation of the CCFR neutrino experiments. Many 
experiments (SLAG, BGDMS, EMG, NMG, E665, GDHSW, BEBG and GGFR) have 
contributed to the heroic and successful effort to obtain a fundamental and precise 
knowledge of properties of partons and of QGD. Gharacteristics of a number of them 



are listed in Tabic 2.1. The fixed-target electron (muon) scattering experiments were 
almost always inclusive, i.e. information on the kinematic variables came only from 
measurements of the incident and scattered leptons. In the charged-current neutrino 
experiments the outgoing muon and the total energy of the produced hadrons are 
measured, and in neutral-current experiments only the latter is detected. 

In deep-inelastic experiments the low x region is correlated with low values of 
Q^, as shown in Fig. 2.1. For fixed-targets the lowest values of x were reached by 



the NMG at GERN and E665 GoUaboration at FNAL applying special experimental 
techniques permitting measurements of muon scattering angles as low as 1 mrad. 
These "small x triggers" and special off-line selection methods were also effective 
against the background of muons scattered elastically from target atomic electrons 
which produce a peak at x =0.000545. Systematic errors on F2 in both experiments 
(in particular those on the ratio of structure functions for different nuclei, -Pf/^l ) were 
greatly reduced as a result of exposing several target materials at the time and/or by 
a frequent exchange of targets in the beam. 



Table 2.1. Fixed-target experiments contributing to the F2 measurements. The x 
and ranges refer to the F2 data; structure function ratio measurements extend 
the low limits by approximately an order of magnitude. 



Beam 


Targets 


Experiment 


Q"(GeV") 


X 


e 


p,d,A 


SLAG 


0.6 - 30 


0.07 - 0.8 


M 


p,d,A 


BGDMS 


7.5 - 230 


0.07 - 0.6 


M 


p,d,A 


NMG 


0.5 - 75 


0.006 - 0.6 


M 


p,d,A 


E665 


0.2 - 75 


0.0008 - 0.6 




Fe 


GGFR, GDHSW 


1.0 - 500 


0.015 -0.6 



In the one-photon-exchange approximation, the differential electroproduction cross 
section is related to the structure function F2{x, Q^) and the ratio R{x, Q"^) of the cross 
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sections for the longitudinally and transversally polarised virtual photons by 



(fa{x,Q'^) Ana 
dQ^dx ^ Q^x 



2 r 



1 



y 



Mxy 
2E 



1 - 



2(1 + i?) 



F2(x,Q2)(2.1) 



where M and m are the mass of the proton and the electron respectively, and E is 
the incident lepton energy. For the HERA kinematics by neglecting M and m this 
expression reduces to: 



d^a 



dxdQ"^ Q^x 



y 



2(1 + i?) 



F2{x,Q^). 



(2.2) 



The function R{x,Q^) has so far been measured only in fixed-target experiments, but 
even here information is scarce. The usual procedure to determine the ^2(2^,^^) is 
to assume a value of R{x,Q^) (theoretical, experimental or a combination of these) 
and then to extract F2{x,Q^) from the data, using an iterative comparison of the 
experimental yield (corrected for acceptance, inefficiency of the apparatus as well as 
for higher-order QED processes) with the electroproduction cross section. For the 
HERA measurements R was calculated using the QCD relation 1^ with the GRV 
structure function parametrization. At small x and the assumed R values can be 
as large as 1. Note that a 20% error on R corresponds to about a 2% uncertainty on 
F2 a.t y — 0.6 for R of about 0.6. The effects due to Z boson exchange for neutral 
current interactions in the presently covered high region for F2 at HERA regime 
amount to a few percent only. New data on R{x,Q^) from fixed-target experiments 
will be discussed in Section 2.5. 

In this Section we present the status of the F2 measurements at HERA and in the 
fixed-target experiments, QCD fits to these data, measurements of R{x, Q^) and tests 
of the (flavour singlet and nonsingiet) sum rules. 



2.2. Measurement of F2{x,Q'^) at HERA 

Both the HI and ZEUS experiments have released new data on structure function 
measurements at small x. It was noted above, that to determine the kinematical 
variables x and Q^, we can use two out of four experimentally accessible quantities: 
the energy, E'^, and angle, 0^, of the scattered electron, and the energy, Eh, and 
average angle, Oh, of the hadron flow. The ultimate method is a global fit of all 
observed quantities, which requires a level of understanding of the detector response 
and of the error correlations that the experiments have not yet achieved. In total four 
methods are currently used in the analyses to reconstruct the event kinematics. The 
electron method, is the method used so far in all fixed-target experiments. Here the 
basic formulae for and y are 

E'^ n 6p „o _ o Op E^ sin'^ 9^ 

= 1 - -f sin^ ^ Ql^ AE'^E, cos^ ^ = ^ (2.3) 

Ee 2 2 I - ye 

The polar angle 6e is defined with respect to the proton beam direction, referred to 
as "forward" region. Ee is the incident electron energy. It remains at HERA the 
most precise way to reconstruct in the whole kinematic range. However at low y 
{y < 0.1) the measurement of x becomes poor and at large y {y > 0.8) the radiative 
corrections which need to be applied to the observed cross section to extract the Born 
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Figure 2.2. Measurement of the structure function with the electron (closed circles) 
and the S method (open squares) by the HI Collaboration. The inner error bar is the 
statistical error. The full error bar represents the statistical and systematic errors 
added in quadrature disregarding the luminosity error (1.5%). 



cross section are very large. The mixed method used by the HI CoUaboration in 1992 
takes from the electron and y fr'om the hadronic variables (yh) according to 

where the sum includes all detected hadrons /i, which have an energy E and 
longitudinal momentum component P^- The resolution of yh is better than the 
resolution of y^ for low y values but becomes inferior at large y values. For the 
double angle method (DA) only the angles of the scattered electron and the 
hadronic system are used. The method is almost independent of energy scales in 
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Figure 2.3. Measurement of the structure function F2{x, Q^) as function of x. Th e 
closed circles are HI data, the open circles are data from NMC U and BCDMS 10|. 
The inner error bar is the statistical error. The full error bar represents the statistical 
and systematic errors added in quadrature disregarding the luminosity errors. The 
curves represent the NLO QCD fit, discussed in Section 2.4. 



the calorimeters but, at very low y, the method is very sensitive to noise in the 
calorimeters. The variables y and are reconstructed from 



yoA 



smt 



;(1 - cos6>h) 



sinOk + sin6e — sin{6e + Oh) 



Qda 



= AEl 



sin6';i(l + cos 9e) 



sin Oh + sin 0^ — sin(0e + Oh) 



with 



Oh 

tan- 



(2.6) 
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This method has been used by the ZEUS Collaboration. A new method used by the 
HI Collaboration |12, 13 , called the S method (E), determines y and from 



where the sum runs over all hadrons in the numerator and over all hadrons plus 
the scattered electron in the denominator. In this method the energy of the incident 
electron in the interaction is reconstructed, which reduces drastically the sensitivity to 
the main radiative process. The resolution in x at low y is good enough to allow the HI 
Collaboration to reach y — 0.01. The resolution at large y is worse but less sensitive to 
radiative corrections than when using only the measurement of the scattered electron. 
For precision measurements of the structure function the different methods are used 



to control the systematics of event smearing and radiative corrections. In Fig. |2.2| the 
results using the electron and S methods are compared for the HI data, showing a 
very good agreement. 

Due to the inevitable beampipe hole for detectors at a collider, the scattered 
electron has to have a minimum angle to leave the beampipe and be detected. From 



(2.3) it foll ows that this leads to a minimum requirement on Q^, which is also visible 
in Fig. 2.1. To study whether Fi still rises at lower values of Q^, several ways to 



increase the acceptance for low were explored. 

The ZEUS detector was improved for electron detection around the beampipe 
by the addition of a scintillator strip detector on the face of the rear calorimeter. 
This allowed the detection of the scattered electron down to smaller angles (175.5° 
compared to 174° in 1993) with a large improvement of the angular resolution (2 
mrad compared to 7 mrad in 1993). The detector was also used to make an event- by- 
event correction to the scattered electron energy arising from the energy loss in the 
inactive material prior to the calorimeter, thereby improving the energy resolution and 
reducing the energy scale uncertainty. Due to the large statistics the electron energy 
scale and angular shift uncertainty for the the HI results have been reduced to 1% 
and 1 mrad respectively. HI could also - due to the excellent accelerator conditions 
at the end of the 1994 running period - make use of the most inner active elements of 
the calorimeter, and increase its low scattering angle acceptance from 173° to 174°. 
In 1993 HI initiated a pilot project to shift the interaction vertex of the collisions 
towards the forward (proton) region. They demonstrated that collisions produced at 
a position shifted by about 70 cm downstream of the detector could still be used for 
Fi analysis, and allowed an increase in the acceptance from 8.5 GeV^ to 4.5 GeV^ 
for the 1993 data. In 1994 about 10 times more data were accumulated using this 
method, and results have been shown by both experiments. 

Another way to access low is by using a sample of deep inelastic events with an 
energetic photon (e.g. E^ >4 GeV) emitted collinear with the incident electron. These 
radiative events can be interpreted as deep inelastic scattering events with a reduced 
incident energy E^ = E,, — E^ which can be reconstructed through the detection of 
the radiated photon in the small angle photon tagger of the luminosity system of the 
experiments. When using the electron method, the "true" kinematic variables yt and 



for such an ep collision are obtained by replacing in (2.3) the nominal beam energy 
by the reduced energy Er. Both experiments have shown data using this process 0|. 

In summary, compared to the previous analyses, the F2 measurement has been 
extended to lower and higher Q"^ (from 4.5 - 1600 GeV^ to 1.5 - 5000 GeV^), and to 
lower and higher x (from 1.8 x 10""^ - 0.13 to 3 x 10"^ - 0.32). The final result is 
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Figure 2.4. Measurement of the structure function F2{x, Q^) as function of for 
the HI, ZEUS, NMC and E665 experiments. The curve is the prediction of the GRV 
parton distributions. 



shown in Fig. 2.3 as a function of x and in Fig. 2.4 as a function of Q^. The error bars 
of the data are now reduced to the 5% to 10% level (except at high and very low Q^). 
The normalization uncertainty has been reduced to 1.5% (2.5%) for HI (ZEUS). The 
rise of F2 with decreasing x is confirmed with higher precision. This rise continues, 
albeit less strongly than at higher Q^, in the region of the lowest available. Scaling 
violations are clearly seen in the plot of F2 versus and will be interpreted in terms 
of QCD below. 

To emphasize the rise of F2 at low Q^, data from the eight lowest bins are shown 
and compared with recent F2 parametrizations in Fig. 2.5. It demonstrates that the 
rise of F2 towards low x is also present in the low region. The measurement is 
in good agreement with recent data from fixed-target experiments E665 and NMC at 
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Figure 2.5. Measurement of the proton structure function F2(x,Q'^) in the low 
region by HI (closed circles: non-radiative events; closed triangles: radiative events), 
together with results from the ZEUS (open squares), E665 (open points) and NMC 
(open triangles) experiments. Different parametrizations for F2 are compared to the 
data. The DOLA and CKMT curves are only shown for the upper row of bins; 
CTEQ3M and MRSG are shown for the lower row; GRV and MRSA' are shown for 
the full range. The values of the ZEUS data shown for the bins = 3.5, 5 
and 6.5 GeV^ are measurements at 3.0, 4.5 and 6 GeV^ respectively. 



higher x values. The curves represent GLAP QCD inspired predictions (GRV fl^ , 
MRS 111 a nd CTEQ and two Regge inspired predictions (DOLA |l7| and 

CKMT [|l8[). The Regge inspired predictions, shown only for the lowest bins, 
are generally below the data. 

The persistent rise of F2 at low x for small indicates that the photoproduction 
regime has not been reached yet. This can also be seen in Fig. 2.6 which shows for 
the low data the strong rise of F2 as a function of W, the invariant mass of the 
7*p system (at low x, W ~ \/Q'^/x). F2 is related to the total cross section of the 
proton-virtual photon interaction (Ttot{l*p) via 

aUl*p)^-^F2{W,Q^). (2.8) 

The F2 growth can be contrasted with the weak rise with W of the total real 
photoproduction cross section in the same range of W, as shown in Fig. 2.(;. The 
different behaviour for = and data at a finite small remains one of the 
interesting questions to be studied at HERA. 
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Figure 2.6. Measurement of the proton structure function F2{W,Q^) as function 
of W^. The inner error bar is the statistical error. The full error represents the 
statistical and systematic errors added in quadrature. Also data , at Q'^ = are 
shown. Note that this figure contains the preliminary HI data ]19| . The curves are 
an ALLM parametrization p^ . 



2.3. Fixed-target structure function data 

The NMC has presented their analysis of the proton and deuteron structure 
functions §], in the range 0.006 < x < 0.6 and 0.5 < < 75 GeV^, as shown 



for F2 in Fig. |2.7| , performed on the almost final sample of events. A clear scaling 
violation pattern with slopes d\nF2/d\nQ^ positive at low x and an "approach to 
scaling" (i.e. a rise of F2 from = to the scaling region) is visible. In this figure 
a comparison of the NMC, SLAC |l| and BCDMS measurements is also shown. 
All three data sets are in good agreement with each other. They were thus used to 
obtain parametrizations of F^ and F2 and their uncertainties, using a 15-parameter 
function. The low x results of the EMC NA2 experiment have earlier been disproved 
by the NMC measurements. The data confirm a characteristic weak x dependence of 
F2 at low Q^, observed for the first time by the EMC NA28 experiment and also 
interpreted in (see also Section 3). 

New measurements of the proton and deuteron structure functions for x >0.0001 



have recently been presented by the E665 Collaboration and are shown in Fig. 2.8 |24|. 
The lowest and x values in their data are 0.2 GeV^, and 8x10"'' respectively. A 
clear pattern emerges from these data at values lower than a few GeV^ , namely a 
weak X, and possibly a stronger than logarithmic dependence, of F2. 
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The E665 region of x, i.e. x >10^^, is now being investigated by both the HI and 
ZEUS Collaborations at HERA. The most dramatic effect visible in the HERA large 

data is a strong increase of F2 with decreasing x, while there is a rather weak x 
dependence of the F2 observed by EMC NA28, NMC, and E665. 

The fixed-target F2 data have had great impact on the determination of parton 
distributions (see e.g. [|5|). It is now seen that these data join well to the results of 
HERA and thus ma ke a joint QCD analysis possible in a large kinematic interval (see 
below). In Fig. ^.9| a detailed comparison of the structure function F2 as function of 
g2 between HI, ZEUS, E665 and NMC is shown for x values around 0.2, 0.07, 0.05 
and 0.012. The data show a smooth continuation over the whole region. It also 
shows a (still) substantially different level of accuracy between the HERA and the 
fixed-target experiments. The former are still expected to improve both in statistics 
and systematics in the next few years. Apart from the above overall agreement, there 
exists however a discrepancy between the NMC and CCFR F2 data at low x (not 
shown). This was much discussed at the Workshop, but did not result in any new 
conclusion. 

Both the NMC and E665 experiments have measured the deuteron to proton 
structure function ratio, F2/F2, extending down to very low values of x. In the case 
of NMC the ratio has been measured directly, i.e. the measurement of the absolute 
structure function is used only for calculation of the radiative corrections. The data are 
usually presented as the ratio F2/F2 where F2 is defined as 2_Pf — F^". This quantity 
would give the structure function of the free nucleon in the absence of nuclear effects in 
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Figure 2.8. Measurements of by the E665 Collaboration. The errors are 
statistical and systematic added in quadrature, a normalisation uncertainty (1.8%) 
is not included. Curves show model calculations of Martin, Stir ling and Roberts, 
Badelek and Kwiecinski, and Donnachie and Landshoff (from |24|). 



the deuteron. The results are presented on the left in Fig. |2.1C p6[ . In both data sets 
the average varies from bin to bin reaching down to (Q / = 0.2 GeV^ at x =0.0008 
for the NMC and (Q^) ^ 0.004 GeV^ at x = 5x10-^ for E665. The results of both 
experiments show that the ratio jF^ remains below unity down to the smallest 
measured values of x. At low x this can be attributed to nuclear shadowing in the 
deuteron predicted to be only weakly x dependent, as observed. It seems unlikely 
that the results can also indicate a difference in F^ of protons and neutrons at low x, 
since e.g. in Regge models the difference between the proton and neutron structure 
functions vanishes with decreasing x ^ . 

New data have appeared on nuclear shadowing. NMC have performed a high 
precision study of the A dependence of nuclear shadowing in the range 0.004 < x < 
0.6 and 1.5 GeV^ < < 60 GeV^. The results are shown on the right in Fig. |2l^ |§. 
These measurements of the ratios F^ jF^ for A = Be, Al, Ca, Fe, Sn and Pb taken 
in conjunction with those on D, He, Li, C and Ca ^ and with earlier data of 
SLAC show a detailed pattern of the x dependence of shadowing. The NMC 
data range from A — 2 to A — 208. The functional dependence of F^ jF^ on A has 
been parametrized as F^ jF^ = cA^"~^'> in each bin of x. The amount of shadowing 
increases strongly with the mass number A. Lower values of x and are covered by 
the nuclear data of E665 in the region x >0.0001 and >0.1 GeV^ Q, as shown in 




Figure 2.9. Detailed comparison of the structure function F2 as function of 
between HI, ZEUS, E665 and NMC, for x values around 0.2, 0.07, 0.05 and 0.012. 
The error bars represent the full errors on the data points. 



Fig. 2.11 (left). A decrease in the amount of shadowing observed by E665 is presently 
under discussion. Shadowing seems to saturate at x about 0.004 as also indicated by 
the NMC data on the F^^Ff and F^ /F^ ratios measured down to x =0.0001 and 
(5^=0.03 GeV^ |p4|. No clear dependence is visible in the E665 data in a wide 



interval of Q , shown in Fig. 2.11 (right), contrary to the preliminary NMC results 
in which positive slopes for the F2'^ / F2 ratio at a; < 0.1 are observed, as shown 
in psf . The shadowing region seems to have another interesting feature: it contains 
a large fraction of large rapidity gap (or diffractive) events, their fraction increasing 
with A Hi. 

2.4- Parton distribution measurements 

Both the HI and ZEUS experiments have performed next-to-leading order (NLO) 
QCD fits based on the Altarelli-Parisi (CLAP) evolution equations on the HERA and 



fixed-target F2 data. Fig. 2.3 shows that the F2 data can be well described by such a 
QCD fit. Note that only data with > 5 GeV^ were used in the fit. The fit result 
was evolved to lower and used as a prediction in the region < 5 GeV^. The low 
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Figure 2.10. Left: E665 results, a) F2 / for three different techniques! 
extracting the ratio. The curve shows a prediction of Badeiek and Kwiecifiski, |27 



b) IFJ^IFl = 2F^/F^ 



1 as a function of x. The NMC data at =4 GeV^ 



are also shown. Errors are statistical The systematic uncertainty is represented by 
the hatched area in Fig. b) (from l26|). Right: NMC results 

on F2 /F^ , averaged over (open symbols) together with earlier results of SLAC 
(closed symbols). Inner/outer error bars represent the statistical/total errors. The 
SLAC-E139 data for silver and gold were used for the comparison with the tin and 
lead data of the NMC, respectively. The photoproduction cross section data are given 
at a small value of x for convenience (from |30|). 



data are found to be well described. This result suggests that, within the present 
accuracy, no (large) higher twist terms are required in this region, contrary to the 
NMC fit, see below. This may need to be reviewed when higher precision data at low 
X become available. 

The scaling violations from the HERA data allow an estimate of the gluon density 
xg{x) at low values of x, while the fixed-target data are needed to constrain the high x 
region. The HI QCD fit [|3 includes only proton data from HI, NMC and BCDMS. 
Additionally the momentum fraction carried by the gluon is imposed to be 0.44. Apart 
from the ZEUS data, the ZEUS fit |Q includes d ata f rom NMC, both on proton and 
deuteron targets. The results are shown in Fig. 2.12 for = 20 GeV^. The error 
bands shown include a careful analysis of the systematics, taking into account the 
correlation between different sources. The results of the two experiments agree very 
well. The resulting gluon distribution shows a clear rise with decreasing x. Similar 
results have been found in ||l^ , which include also other data than those from structure 
functions. In the region x > 10^^ the extracted gluon densities agree with the result 
obtained by the NMC. 

NLO QCD fits have been performed by the NMC to their (earlier) accurate 
measurements of the structure functions and down to low values of x . The 
flavour singlet and non-singlet quark distributions as well as the gluon distribution 
have been parametrized at the reference scale equal to 7 GeV^. All the data with 

> 1 GeV^ were included in the fit. Besides the leading twist contribution a higher 
twist term was also included using a factor 1 + H{x)/Q^ where H(x) was determined 
from the SLAC and the BCDMS measurements [Eol, averaged over the proton and 
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Figure 2.11. Left: E665 results on ja^ . Errors are statistical, the systematic 
ones are marked as shaded bands. Overall normalisation uncertainties have not been 
included. The NMC results have been shown for comparison where available (from 
^). Right: Q2 dependence of the E665 / data (closed 

circles). Errors result from the fit, in which only statistical errors were considered. 
The NMC results are drawn as open squares (from |34[). 



deuteron, and suitably extrapolated to lower values o f x. Results of the QCD fit 
to the proton structure function data are shown in Fig 2.15 and clearly indicate the 
extension of the QCD analysis to the low x and low regions. The contribution of 
higher twists is however substantial at scales of about 1 GeV^. 

Additional information on the gluon density is extracted from hadronic final states. 
The HI experiment has extracted the gluon density from 2-jet events. The method 
and results are given in |4^, |42). The results are limited to the region x > 0.001, and 
agree with the gluon extracted from scaling violations. The E665 experiment has used 
the energy-energy angular pattern of hadrons produced in DIS to extract the gluon 
distribution function of the nucleon, down to a;=0.005. 

An interesting point in the QCD analysis arises from the presence of diffractive-like 
events in the deep inelastic scattering event sample. A contribution of approximately 
10% from events with a large rapidity gap towards the proton remnant, and with the 
characteristics as measured in hadronic diffractive exchange, has been established E^. 
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Figure 2.13. The results of the QCD fit to the F2 data. The sohd hne is the result 
of the QCD fit with higher twist includecLThe dotted curve shows the contribution 
of i^if"^- The errors are statistical (from P9|]). 



These events possibly originate from a different production mechanism than the one 
for the bulk of the deep inelastic data. So far no special account has been taken of this 
in the extraction of partons from inclusive F2 measurements; this topic was discussed 
at the Workshop and should be considered further. 
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2.5. R{x,Q^) measurements 

Extraction of F2{x,Q^) from the data needs information on R{x,Q^). In particular 
the ratio of inelastic cross sections on different nuclei is only equal to the corresponding 
structure function ratio, provided R{x, Q^) is the same for the se nu clei. Results of the 
NMC analysis on R^^-R^ [|| and R'^-RP^, shown in Fig|2]l| (left), demonstrate 
that neither of these quantities exhibit a significant dependence on x and that they 
are both compatible with zero. The NMC reported preliminary measurements of 
R{x, Q^) for proton and deuteron targets as a function of x in the range 0.006 < a; < 
0.14 The average of these measurements ranges from 1.1 GeV^ at the smallest 
X to 15.5 GeV^ at x =0.14. The results show a rise of R with decreasing (small) x. 
Prehminary measurements of R{x, Q^) on a heavy target (Fe), at x >0.01 and > A 
GeV^ (at present) have also been reported by the CCFR neutrino Collaboration p8|| . 
In their data analyses the NMC and E665 experiments assumed R was independent 
of the target atomic mass A and given by the SLAC parametrization |Q valid for 
X >0.1 and > 0.3 GeV^. This parametrization was then extrapolated (with 100% 
error) to 0. Hence there is a need of a theoretical estimate of R (or F^) in the 

region of low x and low Q^. Two ongoing phenomenological studies are expected to 
deliver such estimates soon. In these studies both the perturbative QCD contribution, 
which at low x and low is dominated by the photon-gluon mechanism, and a non- 
perturbative term are taken into account. In the latter contribution is determined 
phenomenologically (Fig. 2.14 (right)) while in Gsj it is fitted to the low data. 



* NMC preliminary 
A SLAC 




Figure 2.14. Left: NMC (preliminary) results R'^ — RP compared with the 
QCD predictions (the curve) and with the results of SLAC (open symbols), from 
|46| . Right: R{x,Q^) in the phenomenological model of 

Badelek, Kwiccinski and Stasto, for two different parton parametrisations (from ]5o|). 

Also at HERA the study of R is an issue of interest. For this purpose it has been 
proposed to operate HERA at lower energies, to have the cross section measurements 
at two y values for a given x, point. Such a measurement could be done in 1996 or 
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1997. A discussion on this topic is presented in Q. 
2.6. Sum Rules 

Several sum rules have been formulated for different combinations of structure 
functions. Strict QCD predictions, valid for — > oo, exist for those involving 
only flavour nonsinglet contributions: the Gross-Llewellyn-Smith and the Bjorken 
sum rules. Experimental measurements of such sum rules provide a stringent test of 
fundamental QCD assumptions. They also in principle permit the extraction of the 
strong coupling constant, as, from the data. Due to the finite of the measurements, 
a predicted value of a sum rule is usually presented in the form of a power series in 
as, the coefficients of which are directly calculated. 

There is no strict QCD prediction for the sum rules containing the flavour singlet 
contributions, i.e. the Gottfried and EUis-Jaffe sum rules. The reason is that singlet 
contributions contain an "intrinsic" dependence. Testing them usually results in 
surprises which teach us a lot about the shortcomings of the simple quark model. 

All the sum rules involve integrations over the whole 0< x <1 interval. This means 
that due to the limited experimental acceptance, interpolations from Xmin to and 
from Xmax to 1 have to be performed. Usually the former is more problematic due to 
a larger contribution of the small x region to integrals and/or to a poor theoretical 
understanding of this kinematic region. Thus the extrapolation a; — > is a major 
source of systematic errors in such sum rule tests. Another source is produced by the 
limited experimental acceptance in at each x value. This usually means that a sum 
rule is measured at a certain Qq, common to all points but at values of Qq which are 
not sufficiently high to exclude a contribution from nonperturbative effects ("higher 
twists" ) . 

2.6.1. Tests of the flavour nonsinglet sum rules. The Bjorken sum rule involves the 
spin structiu'e functions and will be presented in Section 7. The Gross-Llewellyn- 
Smith sum rule, formulated within the Quark Parton Model, states that the integral 
over the valence quark densities is equal to 3 i.e. xF^{x)dx/ x—'i. The QCD 
corrections to this rule have been calculated up to a| The sum rule has been 
tested by the CCFR Collaboration in conjunction with additional low energy data 
from bubble chamber experiments [^^ . The sum rule is fulfilled at the 10% accuracy 
level. The value of as has been obtained at Q^=3 GeV^; it corresponds to as{Mz) 
= 0.108t|]:g!]^(stat.)±0.004(syst.)lg;E]J^^(HT). The uncertainty due to the low scale of 
the measurement (i.e. the presence of "higher twists") dominates the statistical error. 

2.6.2. Tests of the flavour singlet sum rules. The NMC measurements of F2 (fitted 
together with the SLAC and BCDMS data) and of F2 / allow a determination of 
the Gottfried sum i.e. 5*0 = / (F| - F^)dx/x where F| - F^ = 2Fi(l ~ F^/F^)/{1 + 
F^jFl). At Q2 =4 GeV^, neglecting any dependence, Sg was found to be 0.235 
± 0.026 significantly below the simple quark-parton model value of 1/3. This 
is evidence for a flavour asymmetric sea in the nucleon (d sea quarks carry more 
momentum than u), a fact confirmed by the NA51 measurement of the Drell-Yan 
asymmetry in and pn collisions, which gave: u/(i=0.51±0.04±0.05 at a;=0.18 
and Q'^=25 GeV^ Q. Recently a non- negligible Q"^ dependence of — f J* as a 
function of x at low has been reported by the NMC: both the position of the 
maximum and the maximum value of this function change with Q^. This change 
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becomes negligible when higher twist contributions are separated out from the ^Fj 
and F2 / F2 measurements p5[ . The EUis-Jaffe sum rule, involving the spin structure 
functions will be discussed in Section 7. 



3. Low Q^, low X insights from fixed-target data 

3.1. Introduction and basic concepts 

Due to the experimental constraints the fixed-target studies of deep inelastic scattering 
at low X necessarily were correlated with low (Q^ < 1 GeV^). There are two reasons 
why this kinematic region is of special interest. First, as emphasized in section 2, 
the new HERA measurements at small x highlight the importance of a theoretical 
understanding of the connection between the low and high behaviour. The 
second is a practical reason; a unified treatment of low and high is essential for the 
large data analysis, since to implement radiative corrections we require a knowledge 
of structure functions for Q^eas > > 0. 

At low there are constraints on the structure functions Fi{x,Q^) which follow 
from eliminating the kinematical singularities at = from the hadronic tensor 
VF'"'. It is easy to show that as ^ we require 

F2 = 0(Q2) and Fl = 0(Q*). (3.1) 

Hence it is clear that Bjorken scaling, which holds approximately at high Q^, cannot 
be a valid concept at low Q^. 

In dealing with low data we need to introduce the concept of "higher twists" . 
The operator product expansion leads to the representation 

P^M') = E (3-2) 

where the functions Cn{x, Q"^) depend weakly (i.e. logarithmically) on Q^. The various 
terms in this expansion are referred to as leading (n = 0) and higher (n > 1) twists. 
The QCD improved parton model where 

F2ix,Q^) = xJ2 eUq^{x,Q^) + qdx,Q^)] + 0{as{.Q^)), (3.3) 

i 

and which gives approximate Bjorken scaling, retains only the leading twist 
contribution. Physically the higher twist effects arise from the struck parton's 
interaction with target remnants, thus reflecting confinement. For of the order 
of a few GeV^, contributions of the "higher twists" may become significant, see, 
for example, ref. ||39|]. Contrary to the common opinion higher twists, which are 



corrections to the leading (approximately scaling) term (3.3), can only be implemented 



for sufficiently large Q . Thus they cannot correctly describe the low Q (i.e. 



nonperturbative) region since the expansion (3.2) gives a divergent series there. In 



particular the individual terms in this expansion violate the constraint (p.l|). In 



order to correctly describe this region the (formal) expansion (3.2) has to be summed 
beforehand, at large Q^, and then analytically continued to the region of ~ 0. This 
is automatically embodied in certain models like the Vector Meson Dominance (VMD) 
model. To be precise the VMD model together with its generalisation which gives 
(approximate) scaling at large can be represented in a form (3.2) for sufficiently 
large Q^. 
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In practical applications to the analysis of experimental data which extend to the 
moderate values of one often includes the higher twists corrections in the following 
simplified way: 



H(x) 
1 H — 



(3.4) 



where the -F/'"^ is the leading twist contribution to F2 and H{x) is determined from fit 
to the data. This simple minded expression may not be justified theoretically since in 



principle the higher twist terms, i.e. functions Cn{x,Q^) for n > 1 in eq.(3.2) evolve 
differently with than the leading twist term. 

Here we give a brief overview of the parametrizations and the data in the low x, 
low region. We refer the reader to refs. and for a more detailed review of 
the treatment of low problems. 

3.2. Parametrisations of structure functions 

There exist several phenomenological parametrisations (fits) of the structure function 
F2 which incorporate the ^ constraints as well as the Bjorken scaling behaviour 
at large |5^, |l^, 60, ^ |8|. Certain parametrisations |60, ^ |8| also contain 



the (QCD motivated) scaling violations. However, they usually are not linked with 
the conventional QCD evolution. Nor is the low behaviour related to the explicit 
vector meson dominance, known to dominate at low Q^. There exist parmetrisations 
which explicitly contain QCD evolution: |2^, ^ ^ . Most of the parametrisations 
essentially extend the parton model formula for F2 down to the low region 
modifying in a suitable way the parton distributions; the model |^ includes also 
the VMD contribution besides the partonic one and is an absolute prediction (i.e. no 
fitting to the data). The low modifications are typically the following ones: 

(i) Instead of the variable x, modified variables Xi = x{l + Q^i/Q^) are used as 
arguments of the parton distributions where 'i' enumerates the type of the parton. 

(ii) Models which at large include the QCD scaling violations, have the evolution 
in either "frozen" below certain scale Qq which is of the order of 1 GeV^ or 
the evolution in a "shifted" variable + Qq is used. 

(iii) Parton distribution functions are multiplied by form factors of the type 
(5^/((5^ + m|) which ensure vanishing of the structure function F2{x,Q'^) at 
Q' = 0. 

Modifications are absent in the dynamical model which, in principle, is meant 
to describe the structure functions only in the large region even if the QCD 
evolution is extended down to very low scales. Also the recent parametrisation of 
the "parton distributions" uses the variable x instead of x. This model does not 
however extend down to the very low values of (i.e. for < 0.25 GeV^) and, in 
particular, it does not accommodate the photoproduction. 

Parametrisations of F2{x,Q^) differ in their small x behaviour. Most of them 
(except and [^) incorporate at large the steep rise of F2{x, Q^) as the function 
of X with decreasing x which is much stronger than implied (for instance) by the 
expectations based on the "soft" pomeron with intercept ap — 1.08. This steep 
increase of F2{x, Q^) becomes very weak at low Q^. Possible dynamical origin of this 
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effect is different in different models being eitlrer attributed to the absorptive effects 
I p^ , |6^ , |63[ , to the onset of the VMD mechanism or to the pure perturbative QCD 
effects related to the change of the "evolution length" Q. 

There exists practically only one parametrisation of the R{x, Q^ ) structure function 
for the nucleon, i.e. the SLAC parametrisation, based on measurements by SLAC, 
EMC, BCDMS and CDHSW and valid at x >0.1 and >0.3 GeV^. Experimental 
analyses in DIS experiments need to know R down to measured values of x and 
for < < Qmeas- Two phenomenological studies deliver estimates of R in 
the unmeasured region. Both the perturbative QCD contribution, which at low x 
and low is dominated by the photon-gluon mechanism, and a non-perturbative 
term are there taken into account. In the latter contribution is determined 
phenomenologically while in it is fitted to the low data. 



3.3. Experimental data 



The lowest values of x, correlated with lowest values of {x ~10~^ and ~0.001 
GeV^), were reached by the E665 Collaboration at Fermilab by applying a special 
experimental technique which permits the measurement of muon scattering angles as 
low as 1 mrad. At HERA the lowest values of (1.5 GeV^) were recently reached 
by two methods: shifting the interaction point in the proton beam direction in order 
to increase the acceptance of low events and by using radiative events with hard 
photon emission coUinear with the incident electron. Fig. 2. 5, [||. The radiative events 
can be interpreted as non-radiative ones with reduced electron beam energy. 

During the last three years an abundance of new data reaching values smaller 
than 1 GeV^ have appeared. These comprise: results on the proton and deuteron 
structure functions from NMC {x > 0.006, > 0.5 GeV^) @ Fig. and E665 
(x > 0.0001, > 0.2 GeV^) pi Fig. |2!1; results on the deuteron-to-proton structure 



function ratio, F^/F^, measured by NMC and E665 for x > 0.0008, > 0.2 GeV^ 



(NMC |65|1) and x > 0.000 005, > 0.004 GeV^ (E665 ^ 

precise results from these two collaborations on x, A and dependence of nuclear 



), Fig. Qleft), 




Fig. 2.10| (right) and Fig. 2.11; and measurements 
TFig. |2.14|(left), RP and R"^ at low x by NMC [@ and 



shadowing [pO, 
of R^" ~ R^^~l^, i?" - i?P 
i?^ at low X by CCFR Q. 

The above-mentioned data were presented and discussed in Section 2. Here we 
shall add only a few remarks connected with their low behaviour. The data on the 
nucleon F2 display a weak x, and possibly a stronger than logarithmic Q^, dependence, 
at lower than a few GeV^. Observe that also the photoproduction cross section 
between the fixed-target and HERA energies increases rather weakly with energy, cf. 
Fig. 2.6. This should be contrasted with measurements at HERA for larger than 
a few GeV^, which show a strong increase of F2 with decreasing x, Fig. 2.3. The 
QCD analysis of the NMC F2 data, show that the contribution of higher twists 
of the form similar to ( |3.4| ) is moderate even at scales about 1 GeV^ . This is visible 
in Fig. 2.8 where the F^{x, Q"^) at low Q"^ is well described by models ||2^, |l^ directly 
containing higher twist contributions. The higher twists seem also to account for the 
dependence of the F2 — F2 function, cf.section 2.6. 

The F2/F2 ratio, Fig. 2.10| (left), which stays always below unity down to the 



smallest measured values of x, reflects nuclear shadowing in the deuteron, only weakly 
dependent on x. The data are well described by a model , which contains the VMD 
part, essential at low Q^, and which relates shadowing to the diffractively produced 
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final states. The agreement extends over nearly five orders of magnitude in x. No clear 
dependence is visible in the shadowing data in a wide interval of Q^, neither in 
Fi/F^ nor in F^ / F^ |2|, Q, cf. Fig. plll(right), except possibly for F2{Sn)/F2[C) 
at X <0.1, Shadowing thus appears as a leading twist phenomenon. 



34- Outlook 

In this section we have listed the ideas and results concerning the electroproduction 
structure functions in the region of low values of x and . F2{x,Q^) should vanish 
linearly with for ~> (for fixed v), an important property which follows from 
the conservation of the electromagnetic current. The purely partonic description of 
inelastic lepton scattering has thus to break down for low Q^. At moderate the 
higher twist contributions to F2 which vanish as negative powers of are sometimes 
included in the QCD data analysis. One also expects at low that the VMD 
mechanism should play an important role. 

The small x behaviour of F2{x, Q^) is dominated by pomeron exchange. Analysis 
of the structure function in the small x region for both low and moderate values 
of can clarify our understanding of the pomeron. At large the problem is 
linked with the QCD expectations concerning deep inelastic scattering at small x . 
Besides the structure functions (or total cross sections) complementary information 
on the pomeron can also be obtained from the analysis of diffractive processes in the 
electro- and photoproduction. This concerns both inclusive diffraction and diffractive 
production of vector mesons p9| . 

Descriptions of the low Q^, low x behaviour of F2 range from pure fits to 
experimental data to dynamically motivated models. 

There now exists a wealth of measurements of F2 in the low Q^, low x region. 
These include the NMC and E665 results which extend down to very low x and 

and display characteristic "approach to scaling" behaviour, as well as the first 
results from HERA at which extend down to 1.5 GeV^. The data were QCD 
analysed, showing its applicability down to scales of the order of 1 GeV^. Nuclear 
shadowing was studied in great detail for targets ranging from A — 2 to A = 208 
by the NMC and E665 collaborations. Its x, and A dependence were precisely 
measured. Preliminary data on R{x,Q^) have also been reported (NMC and CCFR). 

The fixed-target (unpolarised) structure function measurement programme comes 
to an end in 1996. Many experiments contributed in a great and successful effort 
to learn about properties of partons and strong interactions. Several aspects of this 
knowledge are yet not understood. One of these and perhaps the most challenging one 
is low X dynamics and in particular its dependence on the probing photon virtuality, 
Q^. The new possibilities concerning the study of this problem have opened up with 
the advent of HERA. The data collected there show a very strong increase of F2 with 
decreasing x at high Q^. As the HERA data improve and extend to lower it will 
be informative to see in which region of the strong increase gives way to the slow 



rise evident in photoproduction, see Fig. 2.6. Such data in the transition region open 



up the possibility of a unified understanding of the underlying dynamics. 



4. QCD interpretation 
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4-.1. Introduction 

The behaviour of the proton structure function F2{x,Q^) at small x reflects the 
behaviour of the gluon distribution, since the gluon is by far the dominant parton 
in this regime. Fig. 4.1 shows a sketch of the gluon content of the proton in the 
various kinematic regions. Perturbative QCD does not predict the absolute value of 
the parton distributions, but rather determines how they vary from a given input. For 
instance from given initial distributions at some scale Qq, Altarelli-Parisi (GLAP) Q 
evolution enables us to determine the distributions at higher Q^. GLAP evolution 
resums the leading as \n{Q^/QQ) terms where, in a physical gauge, the \n^{Q/Ql) 
contribution is associated with a space-like chain of n gluon emissions in which the 
successive gluon transverse momenta are strongly ordered along the chain, that is 
< . . . < < Q^. 
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Figure 4.1. The gluon content of the proton as "seen" in various deep inelastic 
{x, Q^) regimes. The critical line, where gluon recombination becomes significant, 
occurs when W ^ 0{as)- W is the ratio of the quadratic recombination term to the 
term linear in the gluon density which occur on the right hand side of the evolution 
equation. 

At sufficiently high electron-proton cm. energy y/s we encounter a second large 
variable, 1/x ^ s/Q"^, and we must resum the leading Q!slog(l/a:) contributions. 
The resummation is accomplished by the BFKL equation In this case we have 
no ordering in q^^^ along the chain, but rather, as we evolve to smaller x, we have 
a diffusion or random walk in Inq^^. The lack of strong ordering means that we 
have to work in terms of the gluon distribution f{x, k"^) unintegrated over the gluon 
transverse momentum kx- As we proceed to smaller x, via the BFKL equation, the 
gluon density / is predicted to increase as x^^ with A ~ 0.5 (on account of the 
increased qxi phase space) and to possess a Gaussian-type distribution in In fc|, which 
broadens as y/\n{l/x). This singular type of growth in x, accompanied by the diffusion 
in Infcp, is the characteristic property of the BFKL gluon density f{x, k^). 

The increase of the gluon density with decreasing x cannot proceed indefinitely. 
Eventually we reach the critical line where gluon recombination effects become 
appreciable, see Fig. 4.1. At the onset these effects can be estimated by perturbative 
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QCD, but finally we enter a region of high density of weakly interacting partons, 
where we have the unusual situation in which we cannot use the normal methods of 
perturbation theory even though as is small — a region of much speculation and 
interest. At higher we can evolve further in x before we reach the critical line, 
since the resolution goes as 1/Q and the transverse area occupied by a parton ^ l/Q^. 



4-2. GLAP-based descriptions at small x 

Traditionally the parton distributions of the proton are determined by global fits to a 
wide range of deep inelastic and related data, see the review in these proceedings |7Q] . 
The starting distributions are parametrized at some scale Qq, typically Qq = 4 GeV^, 
and evolved up in using next-to-leading order GLAP evolution. The parameters 
arc then determined by fitting to all the available data. A wide range of precise data 
exist for x ^ 0.05 and so the partons are well constrained in this region. 

On the other hand for small x, x ^ 10""^, only the structure function F2{ep) is 
measured. In principle these HERA data should determine the small x behaviour of 
the gluon and the sea quark distributions. Roughly speaking the data on F2 constrain 
the sea S and the data on the slope dF2/dlnQ^ determine the gluon g. For example 
if we take 

F^^ xS ^ Asx-^' (4.1) 
9^--9-A,x-^^, (4.2) 

then we might expect to determine A5 and Xg. The most recent HERA data are well 
described by a parametrization with A5 = Ag ^ 0.2 — 0.25 with Qq = 4 GeV^, see 
MRS (A') However, as the data improve w e sh ould be able to determine As and 



Xg independently. To be more precise the slope (O) gives information on Pqg®g. The 
convolution means that the gluon is sampled at higher x. We will see the convolution 
also has implications for the scheme dependence. 

A satisfactory description of the data is also obtained from a set of "dynamical" 
partons evolved up from valence-like distributions at a much lower scale — 0.34 
GeV^ 10. The dynamical model was originally based on the attractive hypothesis 
that at some low scale we just have the u and d valence quark distributions. The 
gluon and sea quark distributions are generated radiatively from the valence quark 
distributions, via the processes q qg and g —^ qq. Unfortunately the radiated 
distributions are spiked towards small x and must be supplemented by inputting a 
contribution at larger x. For instance to describe prompt photon data valence-like 
gluons are needed at the input scale Qq, and to describe the NMC measurements of 
F2 valence-like sea quark distributions are required. It can be argued that Qq = 0.34 
GeV^ is much too low for perturbative QCD to be reliable. The GRV philosophy is 
that GLAP evolution preserves the leading twist and so the parametrization becomes 
physical at some higher scale. The surprise is that the GRV partons appear to give a 
reasonable description as low as « 1 GeV^. 

An idea of the potential steepness of the gluon distribution within GLAP evolution 
can be glimpsed from the double leading logarithm (DLL) approximation which applies 
at small x and large Q^/Qq W^- In this limit GLAP evolution of the gluon gives 



xg{x,Q^) ^ xg(x,Ql) exp 2 



25 1^0 



(4.3) 
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where t = ln(Q^/A^). That is, xg increases faster than any power of ln(l/x), but 
slower than a power of (1/x). Moreover we see that xg increases faster as a; ^ the 
longer is the Qq evolution. The GRV analysis demonstrates that with a flat (or, 

to be precise, a decreasing) input distribution at small x such a rise is obtained with 



the low starting point Qq = 0.34 GeV^. The origin of the growth of (4.3) at small x can 



be traced back to the singular behavior of the gluon anomalous dimension at n = 1 or 



2\ ^, ^-A 



a; = (see below). On the other hand, if the input gluon is singular, xg{x, Qq) ^ x 
with A > 0, as in MRS/CTEQ, then this behaviour overrides the DLL form. In fact 
the present HERA data for F2 cannot really distinguish between these small x gluon 
behaviours. The data are shown in Fig. 2.4 as a function of InQ^ at given values 



of X. From (4^) we see that the slope of these points is a measure of the gluon. 
Also shown on the plot is the description obtained from GRV partons 10 . The GRV 
prediction is a little too steep at small x but this could be easily accommodated by 
a small increase in Qq. In summary the existing data determine only one parameter; 
either A(r:! 0.2 - 0.25) of the MRS/CTEQ global fits or Qg(w 0.5 - 1 GeV^) of the 
GRV-type description. 

In principle, a possible check of whether the rise of F2 at small x is closer to the 



behaviour (4.3) or to a power-like growth may be obtained from the double asymptotic 



scaling behaviour p^. Namely, (4.3) suggests that, instead of the variables ln(l/a;) 



and ln(t/to), one might use the two variables 

a = ^/ln{t/to)Hxo/x), p = v'ln(tAo)/ ln(xo/cc). (4.4) 

For example, by rescaling the observed F2 by a factor R'p which eliminates power-like 



prefactors to the exponential in (4.3) and higher order corrections in the exponent 
the remainder should be proportional to the exponential, i.e. the logarithm of R'pF2 
should be linear in a and independent of p. Data ]7^ show approximate agreement 
with this prediction, provided the scale Qq is chosen to be 1 GeV^ and xq ~ 0.1. The 
resummation oias ln(l/x) contributions will, at sufficiently small x, give effects which 
violate the double scaling form. 

There is, however, some warning against the use of next-to-leading order GLAP 
evolution in combination with a flat input for explaining the observed rise at small x. 
Consider, for simplicity, GLAP evolution for the gluon alone 

dg{x,Q'^) dy f x\ „ 

^59 - diV'Q )• (4-5) 



ainQ2 y ^99 yy 

In moment space this takes the factorized form 

=7..(-,«5)5(-,0^) (4.6) 
where the anomalous dimension in leading order near u> = 

-fgg EE / dXX'^ Pgg{x,as) ^ " (4.7) 

Jo ^ 

In x-space this approximation is equivalent to retaining only the small x approximation 
Pgg w as/x with as = ias/i^- This is the DLL approximation which generates 
a gluon distribution of the form given in ( |4.3| ). The general pattern of the terms 
that occur in the expansion of the anomalous dimensions and the coefficient function 
functions as power series in as and the moment index u) is summarised in Fig. 4.2. 
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Figure 4.2. Possible terms in the perturbative expansion of the anomalous 
dimensions (and coefficient functions). Leading order GLAP and BFKL evolution 
have only the DLL term in common. 



The terms that are included in the full leading (and next-to-leading) order GLAP 
evolution are shown connected by horizontal dotted lines. On the other hand in the 
vicinity of the point to = which governs the small x limit it will be more accurate to 
resum, in the anomalous dimension, all singular terms of the form 

oo 

7.. « E ^« (4-8) 

Tl=l 

which in x-space corresponds to 

n=l 

that is an tj^" behaviour translates into a (logl/a:)""^ behaviour. These terms are 
included in the BFKL equation, and the coefficients An have been computed fz^ . 
Interestingly the coefficients A2 — — — 0. The anomalous dimension jgg has a 
similar structure as jgg, i.e. the leading order consists of terms (as/w)", whereas for 
7gg and jqq the most singular terms are down by one power of a^: as{as/uj)'^. For 
all four cases the coefficients of these leading order terms are known. In particular, 
for the coefficients are non- vanishing (unlike 7gg), positive definite and large fl^ . 
As a consequence the BFKL increasef of F2(~ Pqg g) with decreasing x appears to 
be much more due to the resummation in Pqg (and in the coefficient functions) than 
that for g{'-^ Pgg ® g). It is clear that for a consistent analysis we need higher order 
corrections to 7gg and 7^^: terms of the order as{as/u)Y . It is hoped that they can 
be derived from the next-to-leading order corrections to the BFKL kernel which are 
being computed by Fadin and Lipatov ||77[| . 



H^'he reduction of the BFKL fc^-factorized (4.14) to the GLAP collinear form is discussed in ref. 
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Based upon our present (limited) knowledge of these singular terms of the 
matrix of anomalous dimensions several numerical studies have been undertaken | |78[ | , 
incorporating resummation in various ways. A general conclusion is that GLAP 
evolution with a flat input distribution is rather sensitive to whether resummation 
of the singular terms is included or not. The resummation in 753, although formally 
nonleading, seems to play a key role. On the other hand, if we start from a singular 
input with A ^ 0.3 then the result appears to be rather stable with respect to 



X 



-A 



resummmation. A theoretical issue is the scheme dependence of the resummation 
effects. Catani and Hautmann emphasized at the Workshop that accurate data on 
F]^ could remove the ambiguities and overconstrain the problem. Indeed they showed 
that it is possible to obtain scheme independent evolution equations of the form 



dF2ix,Q^) 
91nQ2 

dFLix,Q^) 



T22 (E>F2+ T2L ® Fl 



Tl2 (E>F2+ Tll ® Fl 



ainQ2 

which inter-relate only physical observables, where the anomalous dimensions are 
computable in perturbative QCD. 

4-3. BFKL-based description of F2 at low x 

This topic is well reviewed in |Q, so we will be brief. The BFKL equation for the 
unintegrated gluon distribution f{x,k^) may be written in the differential form 



df(x, fc|,) ias 2 f '^^T 
dln{l/x) J 



f{x,k!^) f{x,k^) f{x,k'^) 



\k'S - k^l (4/fc^4 + 



(4.10) 



= K®f (4.11) 

where the f{x,k!^) term corresponds to real gluon emissions, and where the 
cancellation of the k'j, ~ k^ singularity between the real and virtual gluon 
contributions is apparent. The gluon distribution /(x, k^) may be calculated by 
integrating down in x from a starting distribution at, say, x — xo — 0.01. The 
BFKL equation effectively sums the leading as log(l/x) contributions; from (4.11) we 
see that as a; — *■ 

/ ^ exp(Alog(l/x)) ~ x-^ (4.12) 

where A is the largest eigenvalue of the BFKL kernel K. For fixed as it can be shown 
lH to be A = (Sas/Tr) 41n2. We need to know the NLO log(l/a:) contributions to 
fully specify the running of as- However, if a physically reasonable prescription for 
the running is assumed then it is found [ [79| that the solution has the form 

/ - C(/c|) x-^. (4.13) 

For running as the integral in (4.10) is weighted more to the infrared (non- 
perturbative) region, but it is found that A ~ 0.5 is less sensitive to the treatment of 
the infrared region than the normalization C. 

The BFKL predictions of the proton structure functions Fi are obtained using the 
kx factorization theorem ||8^, |8^ 



J X J "^T 



(4.14) 
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Figure 4.3. Pictorial representation of the fcy-factorization formula of (4.14). The 
variables x' and kT that are integrated over are, respectively, the longitudinal fraction 
of the proton's momentum and the transverse momentum carried by the gluon which 
dissociates into the qq pair. 



where F^^ are the ofF-shell gluon structure functions which at lowest order are given 
by the quark box (and crossed-box) contributions to photon-gluon fusion, see Fig. 
4.3. In the small x region it is clearly important to tackle the resummation of the 
Q!slog(l/x) contributions and to go beyond (next-to-leading order) GLAP. Indeed 
a satisfactory BFKL prediction of the rise of F2 at small x was obtained prior to 
the HERA measurements fz^ . However, there are ambiguities and limitations in the 
BFKL description of F2. We list them below. 

there is a 



(i) Due to the diffusion of t he gl uon distribution f{x,k^) in 



(ii) 

(iiij 
(iv) 



In k^, 



significant contribution to (4.14) from the infrared region which is beyond the 
scope of perturbative QCD and has to be included phenomenologically. The 
ambiguity is primarily in the overall normalization of F2 rather than the x 
dependence. A physically reasonable choice of the infrared parameters is found 
to give the experimental normalization. 

The BFKL equation only resums the leading order ln(l/a:) terms. The next- 
to-leading order contributions are essential for a stable prediction. Sub-leading 
effects, including energy-momentum conservation on the gluon emissions, are 
expected to suppress the value of A, see, for example, Q. When known, these 
will check whether or not the as(fcy) prescription for the running of as is correct. 
An underlying soft Pomeron contribution has to be included in the small x region, 
determined by extrapolation of the large x values of F2. 

Shadowing corrections to the BFKL equation will eventually suppress the x~^ 
growth. Although they have not yet been fully formulated, the evidence from the 
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observed ratio of difFractive to non-difFractive deep inelastic events, and from the 

persistent rise of F2 at very low w 2 GeV^ , indicates that shadowing effects 

are at most 10% in the HERA regime, 
(v) We need a unified approach which incorporates both the BFKL and GLAP 

resummations. We discuss recent developments below. 
The BFKL gluon distribution / is characterised by a steep rise with decreasing x which 
is accompanied by a diffusion in Ink^. F2 measures only the rise; it is too inclusive an 
observable to probe the kr dependence. For this we need to explore the properties of 
the final state, such as deep inelastic events containing an identified energetic forward 
jet, see section 5. 



4-. 4- Unified evolution 

Because of the possible onset of the BFKL behaviour of F2 in the HERA small x 
domain it is important to study the validity of the GLAP evolution in this region 
and to look for an interpolation between the two approaches. First answers to these 
questions are contained in the studies of resummation effects (see above). Including 
the log(l/a;) resummation effects in GLAP, via the (as/w)" and as{cts /coY^ terms 
in and 7^^ respectively (and also in the coefficient functions) leads to a stronger 
increase of F2 at small x, i.e. we are getting closer to BFKL. 

Ciafaloni, Catani, Fiorani, and Marchesini |84j have proposed a unified evolution 
equation which embodies BFKL evolution at small x and GLAP evolution at larger x 
(see also the review in these proceedings |^). The CCFM equation is based on the 
coherent radiation of gluons, which leads to an angular ordering of gluon emissions. 
Outside the ordered region there is destructive interference between the emissions. For 
simplicity we concentrate on small x. Then the differential probability for emitting a 
gluon of momentum q is of the form 

dP ^ as An- ^0(9-6') (4.15) 
z irq^ 

where successive emissions along the space-like gluon chain occur at larger and 
larger angles. An represents the virtual gluon loop corrections which screen the 



1/z singularity. We can use (4.15) to obtain a recursion relation expressing the 
contribution of the n gluon emission in terms of that of n — 1. On summing we 
find that the gluon distribution satisfies an equation 

F{x,kl,Q^)^F'-°\x,k^,Q^) + 

(4.16) 

J X ^ J ^ 
where F = f /k^. The angular ordering, in the form of the O function, introduces 
a dependence on an additional scale (that turns out to be the hard scale Q of the 
probe) , which is needed to specify the maximum angle of gluon emission. 

When we "unfold" Ar, so that the real and virtual contributions appear on equal 
footing, and then take the leading ln(l/a;) approximation we find that (4.17) reduces 
to the BFKL equation with F independent of (since Q{Q — zq) 1). On the 
other hand, in the large x region Ar ^ 1 and <d{Q — zq) Q{Q — q), which leads to 
GLAP transverse momentum ordering. If we then replace as/ z by Pgg we see that 
(4.17) becomes the integral form of the GLAP equation. 
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Explicit numerical solutions F{x, k^, Q^) of the CCFM equation have recently 
been obtained in the small x region pq |. As anticipated, they show a singular 
behaviour, with A ~ 0.5, a fc-r dependence which broadens as x decreases and a 
suppression of the gluon as compared to BFKL at low . Fig. 4.4 shows a comparison 
of the CCFM fcT-factorization prediction for F2 = F ® F2^ with recent HERA data. 
The DLL dot-dashed curve is obtained by the same procedure except that we replace 
e(Q - zq) by e((3 - q) and set A^, = 1. The difference between the CCFM and DLL 
curves therefore show the value of x at which resummation effects become important. 
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Fig ure 4.4. The continuous curve is the CCFM prediction of F2 obtained in ref. 
]86| compared with preUminary 1994 HERA measurements . The DLL curve is 
oBtained by the same procedure except that we replace S(Q ~ Z(jr\_hy ©(Q — q) and 
sot Aj; = 1. Also shown are the descriptions obtained by GRV ]14[ and MRS(A') 
llSl partons. The figure is taken from ref. |86|. 
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4-5. Open questions 

A completely satisfying explanation of the observed rise at small x cannot avoid 
investigating the transition from the structure function F2 at small x to the 
photoproduction cross section at high energies. Approaching this transition from the 
deep inelastic kinematic region (large and not so small x) where perturbative 
QCD can safely be used, we expect first to see the onset of more complicated 
partonic interactions (e.g. recombination effects) which will screen the strong rise 
of the GLAP or BFKL approximation. So far, these contributions have not been 
observed. Numerical studies based upon the nonlinear Gribov-Levin-Ryskin [^8| 
equation indicate that in the HERA region it will be difficult to identify the screening 
effects: even if they are not negligibly small in the HERA region, they could be 
masked by a change in the input to the GLAP evolution equations. Unfortunately, 
the GLR equations cannot be used to make a reliable quantitative prediction of 
the magnitude of screening effects, due primarily to the lack of knowledge of the 
characteristic transverse area within which the gluons are concentrate d [p9[, but also 
of the size of the enhancement due to QCD multigluon correlations ||90|, |9l|. So at 
present it is not quite clear yet whether recombination effects are relevant at HERA 
and how large they are. 

The conventional GLAP analysis is constrained to leading twist only, and for large 

this approximation is well-justified. However, once we are approaching the low-Q^ 
region, the neglect of higher-twist becomes more and more questionable. Interestingly 
the BFKL Pomeron does include higher- twist contributions, but only a small subset of 
the higher-twist operators that exist in QCD. In the small-x region, among the twist- 
four operators it is the four gluon operator which is expected to become important. 
Evolution equations of higher- twist operators have been studied in | |9^ , and the 
anomalous dimension of the twist-four gluon operator has been calculated in |9^, . 
But so far, the results of these QCD calculations have not been used yet for numerical 
studies. A phenomenological estimate of higher twist contributions at small x and low 

has been done in paj : it indicates that higher- twist effects are surely important 
when less than 1 GeV^, and when increases that they disappear, for large x 
more rapidly than for small x. Further studies of higher twist contributions at small 
X are important and clearly needed. 

It is quite clear that we are still quite far away from a satisfactory understanding 
of the small-x regime. The measurements of F2 at small x have considerably improved 
and indicate a fairly rapid transition from the "soft" x""'"^ behaviour at = to 
a more singular x^°-^ behaviour by ^ 2 GeV^. Acceptable GLAP-based partonic 
descriptions (including next-to-leading order logQ^ terms) are possible down to at 
least ^ 1 GeV^ . Equally well the small x data can be described by BFKL leading 
order a5log(l/x) resummation, via fc-r factorization; though we await the next-to- 
leading log(l/x) contribution for a stable prediction. We must develop a unified 
formalism which embraces both GLAP and BFKL. The CCFM approach is a step in 
this direction. As noted above, more work is needed to quantify the effects of parton 
screening. Of course the partonic formalism must eventually fail as we go to lower 
and lower Q^. It remains a challenge to match this perturbative QCD approach to 
the successful non-perturbative Regge and Vector-Meson-Dominance description of 
proton structure for ~ 0. 
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5. Hadronic Final States 

In order to make further progress to find the appropriate QCD approach to use in 
the small x regime, hadronic final states are analysed. The experiments at HERA 
measure the full hadronic final state, apart from losses in the beampipe. Several 
results on hadronic final states are reported in |Q. In this section three topics will 
be addressed. The first topic is charged particle spectra as measured in the current 
region of the Breit frame. Here one would expect the timclike evolution of the parton 
shower to evolve as in e'^e" annihilation. The second topic is the study of strangeness 
production and a possible discrepancy of the strangeness suppression factor between 
DIS and e+e~ experiments. Strange particle data can also be interpreted in the light 
of possible signals for instanton production. Finally an overview of the hadronic final 
state at low x will be given with a discussion of possibilities of detecting the onset of 
QCD as log 1/x resummation effects. 



5.1. Breit frame: current region 

A natural frame to study the dynamics of the hadronic final state in DIS is the Breit 
frame joj] . In the Quark Parton Model (QPM) the z-component of the momentum of 
the incoming quark is Q/2 before and —Q/2 after the interaction with the exchanged 
virtual boson. In e^e^ annihilation the two quarks are produced with equal and 
opposite momentum, ±y^/2. In the direction of the struck quark ( called the current 
region) in the Breit frame the particle spectra are expected to have no dependence 
on X and a dependence on similar to that observed in e^e~ annihilation jOSt at 
energy s = Q^. 

Both ZEUS and HI |Q have published results of the charged particle scaled 
momentum and multiplicity spectra in the current region of the Breit frame. The 
results indicate these fragmentation properties of the struck quark in DIS are similar to 
those from quarks created in e+e~ annihilation for > 60 GeV^. For < 60 GeV^ 
the mean multiplicity, {rich), tends to be lower in DIS than in e+e^ experiments p8|. 



see Fig. 5.1 



In the low regions measured at HERA it is possible to study the x dependency 
by comparisons with data from previous fixed-target DIS measurements [^9|, which 
are 2 orders of magnitude higher in x than the values probed at HERA. In Fig. 5A it 
is shown that the DIS (rich) data exhibit a marked dependence on x contrary to the 
naive expectation. Also, in this higher x region the agreement between DIS and e~^e~ 
experiments seems to hold. This x dependence can be explained when mass effects are 
taken into account in the DIS kinematics. For (2+1) jet kinematics, the longitudinal 
momentum of the outgoing quark pair is given by p™*^ ^ ^(■^ ~ Q^)/Q where s is 
the square of the invariant mass of the emitted qq or a qg pair. When s, the 

radiation is emitted in the QPM direction, p™'' ^ Pz^^^ ■ However, at low and 
a;, s is likely to be bigger than and the radiation will be emitted in the direction 
opposite to the QPM direction, hence depopulating the current region. Typically the 
smallest resolvable value of s is around 20 GeV^ in the kinematic region of HERA, 
therefore at around 10 GeV^ the emitted radiation has a p™'' = Q/2 or more, 
whereas p^^^^ = —Q/2. 



Mean Charged Multiplicity 
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Figure 5.1. Twice ttie measured ZEUS multiplicity in the current region is compared 
to the results from HRS, MARK I, OPAL, PLUTO and TASSO. Also show is twice 
the current region multiplicity from FNAL DIS PP experiment. All data have been 
corrected for and A decays. 



5.2. Strangeness Production 

The production of strangeness has been studied by the experiments at HERA 



through the production of mesons and A baryons [100, 101 1. These particles 
are reconstructed from the measurements of the charged tracks of the decay particles. 
The acceptance is presently limited to the central region of the detectors, i.e. in 
pseudorapidity rj \r]\ < 1.3, hence the strangeness production measurement mainly 



covers from the current quark fragmentation region. Fig. 5.2 shows K'^ multiplicities, 
corrected for detector effects, for non-diffractive deep inelastic events (i.e. events 
with no visible gap in the direction of the proton remnant, see as function of the 
pseudorapidity rj and the transverse momentum p^. The data are compared with model 



predictions. The MEPS model is based on the LEPTO [102| program and contains 
first order QCD matrix elements, matched with parton showers which approximate 
the higher orders, for the partonic simulation of the hadronic final state. The CDM 
model is based on the ARIADNE [ |103| program and provides an implementation 
of the colour dipole model of a chain of independently radiating dipoles formed by 
emitted gluons (photon gluon fusion events are not described by this picture and are 
added at a rate given by the QCD matrix elements). The calculations are shown 
for standard strange quark suppression of 0.3, i.e. the production of quark-antiquark 
pairs is uu : dd : ss — 1 : 1 : 0.3, tuned essentially on e+e~ hadronic data. It appears 
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that this predicts to many K'^'s, as is seen most prominently in the i] distribution. 
The calculation with a strangeness suppression of 0.2 is in better agreement with the 
data. 

This difference is intriguing and could be a sign of the breakdown of jet universality, 
e.g. due to a "medium dependence" of the fragmentation, or different kind of 



gluon strings in e"^e and ep scattering [104|. However before speculating on such 



possibilities, additional data are needed. It should be checked whether it is the K° 
rate which is lower, or whether all particle rates are affected, e.g. by studying K/tt 
ratios. Other strange particle species, especially anti-baryons, should be studied as 
well. 

Strange particle production is also of interest for the study of QCD instantons. 
Instantons are related to non-perturbative effects of QCD, a nd e ffects on the hadronic 
final state in ep collisions at HERA were recently calculated |105|. Due to the expected 
isotropic decay of a dense partonic system, events with an unusual structure are 
expected. The uncertainty on the expected instanton rate is however large. Due to 
the flavour symmetry in the instanton decay, less strangeness suppression is expected, 
leading to an average density of 0.4 K^s per unit of pseudorapidity in the region of 
the instanton [106|. The shape of the rapidity distribution is expected to be different 
from DIS events without an instanton. Clearly enhanced strangeness production will 
be one of the features to search for instantons, others being large Et in a limited 
rapidity interval, isotropically spread energy, etc. Searching for instantons is one of 
the exciting topics of hadronic final state studies at HERA. 



5.3. Tests of QCD in the small x regime 

In Section 2 it was shown that the the proton structure function exhibits a strong 
rise towards small Bjorken-x. Originally, this rise caused much debate on whether the 
HERA data are still in a regime where the QCD evolution of the parton densities can 
be described by the CLAP Q evolution equations, or whether they extend into a new 
regime where the QCD dynamics is described by the BFKL evolution equation, 
but it was pointed out in Section 4 that present Fi measurements do not yet allow us 
to distinguish between BFKL and conventional CLAP dynamics, and are perhaps too 
inclusive a measure to be a sensitive discriminator. Hadronic final states are expected 
to give additional information and could be more sensitive to the parton evolution. 

For events at low x, hadron production in the region between the current jet and 
the proton remnant is expected to be sensitive to the effects of the BFKL or CLAP 
dynamics. At lowest order the BFKL and CLAP evolution equations effectively resum 
the leading logarithmic aslnl/x or aglnQ^ contributions respectively. In an axial 



gauge this amounts to a resummation of ladder diagrams of the type shown in Fig. 5.3. 
This shows that before a quark is struck by the virtual photon, a cascade of partons 
may be emitted. The fraction of the proton momentum carried by the emitted partons, 
Xi, and their transverse momenta, /cti, are indicated in the figure. In the leading log 
CLAP scheme this parton cascade follows a strong ordering in transverse momentum 
^Tn ^ ^Tn-i ^ •■• ^ ^Ti, whilc there is only a soft (kinematical) ordering for the 
fractional momentum Xn < Xn-i < •■• < a^i. In the BFKL scheme the cascade follows 
a strong ordering in fractional momentum Xn <C a^n-i ^ ••■ ^ a^i, while there is 
no ordering in transverse momentum. The transverse momentum follows a kind of 
random walk in logfcy space: the value of kTi is close to that of fcTi-i, but it can be 
both larger or smaller. 
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Figure 5.2. Corrected multiplicities in non-diffractive deep inelastic scattering 
events as a function of rj (top) and (bottom) compared to the predictions of the 
MEPS and CDM Monte Carlo program. The inner error bars are statistical and the 
outer error bars contain the full statistical and systematical errors. 



Several measurements of the hadronic final state have been suggested to exploit 
this difference at the parton level. The idea is to find observables which may reflect 
both of the BFKL characteristics of the unintegrated gluon distribution f(x,k'^), that 
is the x~'^ growth and the diffusion in log fc^, as x decreases. Here we discuss those two 
measurements for which new information was provided during the Workshop: first, 
in section 5.4, we study the transverse energy (Et) flow in the region away from the 
current quark jet and second, in section 5.5, we investigate the distribution of deep 
inelastic events containing an identified forward jet, that is, a measured jet as close as 
possible to, but distinct from, the proton remnants. 

Apart from numerical calculations, predictions for final state observables are also 
available as Monte Carlo models, based upon QCD phenomenology. The CDM model 
description of gluon emission is similar to that of the BFKL evolution, because the 
gluons emitted by the dipoles do not obey strong ordering in kr |107|. The CDM does 
not explicitly make use of the BFKL evolution equation, however. The MEPS model 
is based on CLAP dynamics; the emitted partons generated by the leading log parton 
showers are strongly ordered in kx- 



5.4- Energy flow in the central region 

Due to the absence of kr ordering the BFKL approach is expected to give a larger 
transverse energy, Et , in the hadronic final state in the central region of the hadronic 
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Figure 5.3. Parton evolution in the ladder approximation. The selection of DIS 
events containing a forward jet is illustrated. 



centre-of-mass than the GLAP approach at low x. This corresponds to the very 
forward region of the detectors in the HERA laboratory frame. In this central region of 
the hadronic center-of-mass Golec-Biernat aZ , |108] show using the GLAP approach 
that the partonic ), increases with increasing x, while for BFKL the 

reverse is true, {Et) decreases with increasing x. It should be noted though that the 
Et spectra are difficult to calculate in the GLAP framework because the Et weighting 
emphasi zes unsafe regions of phase space. 

Fig. bA shows the transverse energy flo w as a function of pseu dorapidity, -q, in the 
laboratory frame as measured by both HI [109| and ZEUS |11C|. The level of Et is 
almost flat at « 2 GeV/(unit of rf). Also shown are the partonic calculations from 
which seem to indicate a preference towards the BFKL approach over GLAP. 



Fig. shows the mean transverse energy in the pseudorapidity range (as measured 
in the hadronic center-of-mass system) —0.5 < rj* < 0.5 by HI [ |109[ . The data 
exhibit a rise in (Et) with decreasing x as can also be seen from the BFKL partonic 
calculations shown in the same figure. To definitively comment on whether these plots 
are signatures for BFKL it is necessary to understand the underlying hadronisation. 

In order to investigate hadronisation effects in this region two Monte Carlo 
generators have been studied: the CDM model (labeled ARIADNE 4.07 on the figure) 
and LEPTO 6.1, both of which are based on the Lund string fragmentation framework. 
The data shown in Figs. |3.4| and |3.5| are reasonably described by CDM, whereas the 
overall Et predicted by LEPTO 6.1 is far too low. Unfortunately this version of 
LEPTO is very sensitive to the cut-off applied to avoid divergences in the matrix 
element — the boson-gluon-fusion (BGF) process, which is the dominant 0{as) graph 
at low and x, creates two Lund strings as opposed to only one from the QCD 
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Figure 5.4. Mean transverse energy per unit rapidity in the forward region of the 
HERA lab. system. 



Compton process or a ©(a^) scattered quark. The effect of the presence of strings 
is an increased amount of Et available in the hadronisation phase. A newer version 
of the program, LEPTO 6.3, reduced this dependence on the matrix element cut-off 
by treating the scattering off a sea quark in a manner similar to the 2 string scenario 
of the BGF case. Using this version of the generator allows a better description of 
the data, though there are still some problems in describing the x dependence in the 
lowest region as shown in Figs. 5.4 and 5.5. However, the measurement seems to 
be rather sensitive to the non-perturbative hadronisation phase and can, in its present 
form, no longer be considered as a direct sensitive test of BFKL dynamics. 



5. 5. Deep inelastic events containing a forward jet 



Another possible signature of the BFKL dynamics is the behaviour of deep inelastic 



{x,Q ) events which contain a measured jet {xj,kj^j) in the kinematic regime where 
5^ (so as to neutralise the ordinary gluon radiation which would have arisen 



/r2 



from GLAP evolution) and where the jet has longitudinal momentum fraction xj as 
large as is experimental feasible (xj ~ 0.1). The aim is to observe events with x/xj 
as small as possible. According to BFKL dynamics the differential structure function 



has a leading small x/xj behaviour of the form |111 



9 + g (9 + 9) 
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Figure 5.5. Transverse energy as a function of x for three different values of Q^. 
Tlie transverse energy is measured in the CMS in the pseudorapidity range —0.5 to 
+0.5. 



where the parton distributions are to be evaluated at {xj , k^j ) — where they are well 
known from global analyses. The idea is to see if the DIS + forward jet data show 
evidence of the {x/xj)~^ behaviour. Jets are generally expected to be more robust 
against hadronisation effects than is the Et flow. 

We have studied DIS events at small x which have a jet with large xj. A cone 
algorithm is used to find jets, requiring an Et larger than 5 GeV in a cone of radius 
R = a/At^TA^ = 1.0 in the space of pseudo-rapidity rj and azimuthal angle 4> in 
the HERA frame of reference. In order not to confuse the forward jet with the one at 
the top of the ladder the requirement y > 0.1 was imposed to ensure that the jet of 
the struck quark is well within the central region of the detector and is expected to 
have a jet angle larger than 60°. Experimentally a cone algorithm is used to find jets, 
requiring an Ex larger than 5 GeV in a cone of radius R = \J A?/^ + Ac/)^ — 1.0 in 
the space of pseudo-rapidity ry and azimuthal angle </> in the HERA frame of reference. 
Jets are accepted as forward jets if x^ > 0.025, 0.5 < k^j/Q^ < 4, 6" < 9j < 20° and 
kxj > 5 GeV, where 9j is the forward jet angle and kxj is the transverse momentum 
of the jet. 

These selection criteria allow a study the cross section of forward jet production 
in the region w 20 GeV^ and 2- lO"** < a: < 2- 10"'^. Hence the ratio xj/x is always 
larger than 10. 

The resulting number of events observed with at least one forward jet in the 
kinematical region 160° < 9e < 173° and E^ > 12 GeV is given in Table 5.1 and 
compared to expectations of the MEPS and CDM models after detector simulation 
and corrected for background. The measured cross section for forward jets satisfying 
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X range 


data 
events 


MEPS 
events 


CDM 

events 


cr(ep jet+X) 
(Pb) 


2 • 10-'' - 1 • 10-'' 


271 


141 


282 


709 ± 42 ± 166 


1 • 10-^ - 2 • 10-^ 


158 


101 


108 


475 ±39 ±110 



Table 5.1. Numbers of observed DIS events with a selected forward jet, corrected 
for radiative events faking this signature. These may be directly compared with 
the expectations from the Monte Carlo models. The measured cross section ep — ► 
jet+X for forward jets is also given. The errors reflect the statistical and systematic 
uncertainties. 



the cuts given above is also presented in Table 5J. It has been corrected for detector 
effects using the CDM. The ratio of the jet cross section for the low x to the high x 
bin is 1.49 ±0.25. 
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Figure 5.6. Comparison of the BFKL calculation and the approximate analytical 
calculation of the three-parton matrix elements, compared with the data of the HI 
experiment. 



In Fig. |5.6| the data are compared with a recent calculation |112 . The calculation 
(at the parton level) has exactly the same cuts as the measurement. The solid line is 
a BKFL calculation, while the dashed line is an approximate analytical three-parton 
matrix element calculation without the BFKL ladder. The results of HI are compared 
with the calculations. The BFKL curve agrees well with the data, while the lower 
curve is in clear disagreement. Several corrections to this result should be considered, 
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as discussed in [112] ( i.e. the calculation is at the parton level while the data are 
at the hadron level). Nevertheless the agreement with data is encouraging and such 
measurements should be explored further with larger statistics event samples as they 
have a large potential to reveal BFKL effects in the data. 

A final general comment is in order. The smaller the value of xjxj the larger 
is the BFKL effect and the more dominant is the leading log(l/a;) formalism. As 
with all BFKL predictions, the reliability can only be quantified when the sub- 
leading corrections are known. Moreover by measuring properties of the final state we 
inevitably reduce the "reach" of HERA. For example in the present case we require 
xjxj to be as small as possible, yet experimentally jet recognition demands x^ ^ 0.1, 
and so we lose an order of magnitude in "reaching small a;" . 



6. Charm production 

As we have seen above the gluon, the dominant parton at small x, has the least well 
determined distribution. Due to the theoretical uncertainties it is difficult to make 
unambiguous determinations of the gluon from the scaling violations that are observed 
in the HERA measurements of F-i. The longitudinal structure function Fj^ gives, in 
principle, a much more direct determination of the gluon. In practice a sufficiently 
accurate measurement of -Fl, which requires different beam energies, is exceedingly 
difficult. Here we note that F2 itself contains an appreciable part which is directly 
sensitive to the gluon density, namely the contribution arising from charm production 

F2- 




Figure 6.1. Predictions for compared with preliminary HI [113 data which 
became avai lahle after t he W orkshop. The curves show the predictions from MRS(A') 
partons |15| and from W\M using GRV partons (described in section 6.1), together 
with the CCFM prediction for rric = 1.4 and 1.7 of rcf. [ p6| . The figure is adapted 
from ref. feal. 



At this Workshop both Vogt and Sutton presented predictions for F^. The 
predictions are summarized in Fig. 6.1, which also contains preliminary measurements 
of F2 that became available after the Workshop | |113{ . The figure shows two 
calculations based on GLAP-evolved partons, (GRV, MRS), together with a prediction 
obtained using the unintegrated gluon distribution obtained by solving the CCFM 
equation (see section 4.4). All three calculations are quite different. MRS treat charm 
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as a parton and set c{x, Q^) = below some threshold < rn? and then use NLO 
GLAP evolution with rric = 0. The calculation |114] using GRV partons treats charm 



as a heavy quark, not as a parton, and uses fixed-order (NLO) perturbation theory. We 
discuss this approach in detail below. Finally the CCFM predictions are obtained 
from the unintegrated gluon distribution using the fcT-factorization theorem. They 
sum the LO log(l/a;) contributions, but with the angular ordering constraints imposed. 



6.1. Fixed- order prediction for F2 

At this Workshop, Riemersma presented the NLO QCD formalism for the cc 
contribution to F2. It is given by 115 
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■2,3 



- 47ras(^^) 
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X y 



'2,1 



mi 



(6.1) 



where here the sum is over the light quark flavours (rif ~ 3). The (MS) mass 
factorization scale fi has been put equal to the renormalization scale entering the 
strong coupling as. The lower limit of integration over the fractional initial-parton 
momentum y is given by j/min — ax = {1 + 4r7i^/Q^)a:, corresponding to the threshold 
s — 4m^ of the partonic ce nter- of-mass (cm.) energy squared. The coefficient 
functions C2 can be expressed [115] in terms of the dimensionless variables 
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4to? 



-i=|r^-ii-i 



(6.2) 



In leading order, 0{as), F2 is directly sensitive only to the gluon density g{y, /i^) 
via the well-known Bethe-Heitler process 7*5 ^ cc [116 . A comparison of the 



various contributions to F2 in NLO shows that for the physically reasonable scales /z, 
^ ~ 2mc, . . . \JQ'^ + 4rn2 (gee below), the quark contribution in (6.1) — which is not 
necessarily positive due to mass factorization — is very small, typically about 5% or 
less. Therefore F^ does represent a clean gluonic observable also in NLO. 

Before considering the phenomenological consequences of ( 6.1), it is instructive 
to recall the 77-dependence of the gluonic coefficient functions |ll5| , 117 1. They are 
displayed in Fig. 6.2 for two values of typical for deep-inelastic small- a; studies 
at HERA, Q2 = iq, 100 GeV^. Here and in the following we take = 1.5 GeV. 
The comparison of the NLO coefficients and c*^^-* with the Bethe-Heitler result 
c(o) reveals that potentially large corrections arise from regions where c^°^ is small, 
namely from very small and very large partonic cm. energies. These corrections 
are due to initial-state-gluon bremsstrahlung and the Coulomb singularity at small s 
(small 77), and due to the flavour excitation process at s » 4m^ (?7 » 1). The large 
logarithms ln[s/(4m^)] originating from flavour excitation have been resummed — at 
the expense of losing the full small-s information of (6.1) — by introducing a charm 
parton density, leading to the so-called variable-flavour scheme [118|. The importance 
of these corrections in the HERA small- a; regime considered here will be investigated 
below. 

The first question to be addressed in order to judge the phenomenological 
usefulness of F^ as a gluon constraint is whether or not the available NLO expression 
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(6.1) is sufRcient for obtaining results which are stable under variation of the 
(unphysical) mass factorization scale. It has been argued in |119| that one should 
use /X ~ 2mc- This is motivated by the conjecture that /i is controlled by s and the 
fact that the integrand in (6.1) is maximal close to the lower limit, s ~ 4m^. The range 
of significant contributions in s, however, broadens considerably with increasing Q^, 
as we shall see below, hence fj, ~ \/Q'^ + 4m^, chosen in |115, |ll(| , appears at least 
equally reasonable. Therefore, we will estimate the theoretical uncertainty of F2 in 
NLO by varying the scale between /i = nic and fj, = 2yJ Q"^ + 4m^. The corresponding 
results for are shown in Fig. 6.3 for some fixed values of x and Q^, using the parton 
densities of GRV jlj], which are appropriate here as they treat charm as a massive 
quark and not as a parton. At small-x, x < 10~^, the scale variation amounts to at 
most ±10%. Also displayed in Fig. 6.3 are the results for the three CTEQ2 parton 




Figure 6.2. The rj-dependence of the coefficient functions as parametrized in [Ll£| 

^2,9 

given in the MS scheme. 



for two values of with rtic = 1.5 GeV. The scheme-dependent quantity ci „' is 



sets |]120|| at two selected values of x. Obviously the scale stability at small x does 
not significantly depend on the steepness of the gluon distribution. Consequently, the 



NLO results of |115, 116 seem to provide rather sound a theoretical foundation for a 
small- a; gluon determination at HERA, despite the large total cm. energy s ^ 4m^ 
which might suggest a destabilizing importance of ln[s/(4m^)] terms. At large x, 
X « 0.1, on the other hand, the scale dependence of F2 is rather strong, especially 
at low Q^. Here, quantitatively reliable results require further theoretical input by 
resumming the large small-ry threshold contributions mentioned above. However, F2 
is small in this region. 

The next issue we investigate is the locality (in the momentum fraction y) of the 
gluon determination via Fg'^. For this purpose, we replace the upper limit in (6.1) by 
a varying maximal gluon momentum j/max < 1- The contribution from initial-parton 
momenta smaller than ymax to F^fix^Q^), denoted by F2 {x , Q'^ , Umax) , is presented 
in Fig. 6.4 for the GRV parton distributions |l4j. At scales /z ~ ^Q'^ + 4m^, about 
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80% of F2 originates in the region y,nin — o,x < y ^ iUmin- 
very similar for the CTEQ2 parton densities |12C|. Thus F2 



Again the situation is 
aUows for rather focal 



a determination of g{y,fj,^ ~ x/Q^ + 4to^). The partonic cm. energies in the region 
contributing 80% to the structure function F2 are given by 77 ^ 5 (20), corresponding 
to s < 60 (180) GeV^, at ^ iq (iqo) GeV^, respectively. This imphes that for the 

values under consideration here, the plateau region of and C2 at large 77 (c.f. 
Fig. 6.2) does not play any important role. Resummation of large- 77 flavour excitation 
logarithms is thus neither necessary nor appropriate, especially if it implies additional 
approximations in the more important small-s region |118|. A similar observation has 
already been made in 119 for /i ~ 2mc at low Q^. Note, however, that the latter 
scale choice leads to a considerably wider important range of s at high Q^, somewhat 
in contrast to its motivation described above. 

The expected x-dependence of F2 and its relative contribution to the total proton 
structure function F2 are displayed in Fig. 6.5. F2 is large in the HERA small- a; 
region (see Fig. 6.1). The absolute magnitude amounts to about 0.2 . . . 0.4, i.e. F2 is 
as large here as the total F2 in the BCDMS/NMC fixed-target kinematical regime, 
making up up to a quarter of F2 measured at HERA. This is in contrast to the bottom 
contribution F2, which reaches at most 2 . . . 3%. The size of F2 renders a reliable 
(fully massive) treatment mandatory in any precise analysis of F2 at small x. The 
sensitivity of F2 to the gluon density is illustrated by the difference of the CTEQ 2MF 
(flat xg{x,^J? = 2.6 GeV^) and the GRV (steep gluon) ||l|l expectations. There 

is quite some discriminative power of F2 , especially close to the lower end of the 
range considered here, ~ 10 GeV^. Moreover, by measuring up to about 100 GeV^ 
in , the rapid growth of yg{y ^ 1,A*^) predicted by the GLAP equations can be 
rather directly tested down to y ~ 10~^. A theoretical obstacle to an easy accurate 
gluon determination via the charm structure function is the dependence of F2 on the 
unknown precise value of the charm quark mass mc- A 10% variation of rric to 1.35 
GeV and 1.65 GeV as also considered in Fig. 4 leads to a ±15 . . . 25% (5 . . . 10%) effect 
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Figure 6.3. The dependence of on the mass factorization scale ^ in the region 
TOc < u. < 2(Q'^ -I- 4m^)'''^ at sel ecte d values of x and . The parton densities of 
GRV |l^(at all x) and CTEQ2 [[2^ (at x = 10-^,0.04) have been employed. 
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Figure 6.4. The contribution of the initial-parton momentum region ax < y < y' 



Ft' 

to F2 at small-x for two choices of the scale fi, using the parton densities of ]14[ . The 
arrows indicate the values of y™''" at which 80% of the complete results are reached 
for fj. = (Q2 +4m2)0-5. 



at = 10 (100) GeV^, respectively, with respect to the central curves. 
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Figure 6.5. The x-dependence of and /F2 at some fixed values of Q^, as 
expected from the GRV gl ijnn density jlH . Also shown are as obtained from 
the CTEQ 2MF parton set ]120| and the charm mass dependence of the predictions. 
fj, = (Q^ + Am'^y^-^' has been employed. 
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6.2. Outlook for 

The measurement of F| at HERA should serve as a sensitive probe of the gluon at 
small X. In section 6.1 we have seen that the fixed-order QCD predictions are stable, 
with scale variations of less than ±10%, and offer a rather local measurement of the 
gluon. Flavour excitation contributions from s 3> 4m^ become important only at 
scales higher than those relevant for the small x observations at HERA. 

Investigations of the effects of log(l/a;) resummations have started (see Fig. 6.1 
and also |6|) but as yet the theoretical framework is much more incomplete than 
the fixed-order approach. Ultimately these effects may have to be incorporated in a 
detailed analysis of future precision data for 



6.3. J /ip production as a probe of the gluon 

It has long been advocated that inelastic J/ij] photoproduction at HERA may serve 



as a measure of the gluon — see, for example, [121 or the review |12S] which consider 
the process at leading-order accuracy. Recently the higher-order corrections to this 
process have been calculated |123, 124 1. A detailed analysis of the spectra in the high 
energy range at HERA shows that the perturbative calculation is not well-behaved 
in the limit pt — > 0, where px is the transverse momentum of the J /ip- No reliable 




Figure 6.6. Total cross section for inelastic J/4> photoproduction as a function of 
the photon-proton energy for different^wametrizations of the part on distribution in 



the proton. Experimental data from [125|. The figure is from |124 



prediction can be made in this singular boundary region without resummation of large 
logarithmic corrections caused by multiple gluon emission. If the small px region is 
excluded from the analysis, the next-to- leading order result accounts for the energy 



dependence of the cross section and for the overall normalization, see Fig. 6.6 |124|. 
However, since the average momentum fraction of the partons is shifted to larger 
values when excluding the small-p_L region, the sensitivity of the prediction to the 
small-x behaviour of the gluon distribution is not very distinctive. 

Diffractive J/ip photoproduction appears to offer a more promising way to 
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Figure 6.7. The measurements of the cross section for diffractive J/i/i 
photoproduction compared with the full perturbative QC F) pr ediction obtained from 
the three latest sets of partons. The figure is taken from |128!. 



distinguish between the gluon distributions. Since this is essentially an elastic process 
the cross section is a me asure of th e square of the gluon density. To leading order the 



cross section is given by 126 127 
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^2, 



asiQ ) 



(6.3) 



X — M^/W'^, where W is the 7p cm. energy. In a study which 



with Q = iAf| and 

originated at the Workshop |125], corrections to this formula have been calculated 
and comparisons with HERA data made, see Fig. 6.7. It was emphasized that the 
W dependence, rather than the normalisation, was the more reliable discriminator 
between the gluons. The power of the method is evident from Fig. 6.7, which appears 
to favour the MRS (A') gluon. 



7. Spin physics 

Interest in spin phenomena in deep inelastic scattering revived in the eighties after the 
European Muon Collaboration discovered that the spin dependent structure function 
of the proton violates the Ellis- Jaffe sum rule and that the quarks probably carry only 
a small part of the total proton spin. The problem of the origin of the proton spin has 
led to an intense experimental and theoretical activity but in spite of that it has not 
yet been answered conclusively. Here we briefly review the status of the spin effects 
in deep inelastic scattering; experimental results and their theoretical interpretations. 
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Finally its future prospects will be outlined. 



7.1. Cross section asymmetries and sum rules 

The deep inelastic lepton-nucleon scattering cross section is the sum of a spin 
independent term a and a term proportional to the lepton helicity hi: 



(7.1) 



Only longitudinally polarised Icptons will be considered and the spin vector s; is 
thus related to the lepton four momentum vector k. Act gives only a small contribution 
to the total deep inelastic cross section. It depends on the two structure functions g\ 
and fif2 and can be expressed as 



Act = cos^) Actii + sin'^ cos4> Acr_L, 



(7.2) 



with 



d^Aaii 
da;dQ2 

d-'^AcjT 
d^dQ2d^ 



167ra^j/ 



^2y2 



-92 



= — COS 




(7.3) 



In the above expressions, ip denotes the angle between the lepton and the nuclcon 
spin and </> the angle between the scattering plane and the spin plane; furthermore 
Aaj_ = Aar/ coscj). y = v/E and 7 = IMxj \fQ^ are the usual kinematical factors. 
In experimental measurements, two asymmetries can be defined: 



A 



Aaii 



and 



Acri 



(7.4) 



2ct 2ct 
These asymmetries are directly related to the virtual photon asymmetries, A\ and Ai: 



A|| = D{Ax + 77^2), A^ = D{A2 - ^Ai), (7.5) 

where 

. 91- 1^92 . 91+92 

Ai = — , A2 = 7— pT— • (7-6) 

D, often called the dcpolarisation factor of the virtual photon, depends on y and the 
structure function R = Fl/Ft; factors rj and ^ depend only on kinematic variables. 
Ai and A2 are often interpreted as virtual photon - nucleon asymmetries. They satisfy 
the boimds \Ai\ < 1, IA2I < VR. 

Within the QPM the spin dependent structure function gi is given by 

9i{x) = lY.e^[Aqi{x) + Aqi{x)], (7.7) 

with Aqi (x) = qf (.x) — q~ (.x) , where q^ are the distribution functions of quarks with 
spin parallel (antiparallel) to the nucleon spin. Less obvious is the meaning of 52 
which contains a leading twist part, completely determined by g\ and a higher twist 
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part, the meaning of which is subject to debate [129|. In QCD, gi evolves according to 
Altarelh-Parisi equations, similar to the unpolarised ones. Corresponding coefficient- 
and splitting functions have recently been calculated up to order a| [130|, permitting 
the next-to-leading order QCD analysis of gi and thus a determination of the polarised 
parton distributions, Aqi(x,Q^). Various groups have used the recent data to 
determine these distributions, taking into account the leading order QCD correctionsf . 
A comparison of different parametrizations | |132| shows that the polarised valence 
quark distributions Aut,(x, Q^) and Adv{x, Q^) can be determined with some accuracy 
from the data, while the polarised sea quark and gluon distributions Aq{x, Q^) and 
Ag{x, Q"^) are only loosely constrained by the structure function measurements. 

Contrary to 171 and 52 > definite theoretical predictions exist for the first moment of 
(71, Fi = gi{x) Ax, which measures the expectation value of the axial vector current 
between two nucleon states. Two sum rules exist for Fi. The fundamental one was 



obtained by Bjorken |133| from the current algebra and isospin symmetry between the 
proton and the neutron: 



9A_ 

9v 



{Au - Ad) 



(7.8) 



where gA and gy are the axial and vector weak coupling constants in the neutron beta 
decay. The QCD corrections to this sum r ule h ave been computed up to the order 



134 1 and the 0{ag) have been estimated 135 



Separate sum rules, obtained by Ellis and Jaffe |136(] , hold for the proton and the 
neutron: 



F- 



,p(n) _ 



1 

± — 

12 



gA 
gv 



36 ^ 9 



(7.9) 



Here AS = Au + Ad + As is the flavour singlet axial coupling and Aq denote first 
moments of the spin dependent parton distributions in the proton, Aq = Aqi(x) dx\ 
as (and \gA/gv\) are related to the symmetric and antisymmetric weak flavour-SU(3) 
couplings in the baryon octet. If the flavour-SU(3) is exact then og can be predicted 
from measurements of hyperon dacays. There is however no prediction for AS, except 
for As=0. In this case AS = ag, as was assumed in the original formulation by Ellis 
and Jaffe Jl36|. QCD corrections to these sum rules have been calculated up to the 
order a| |137| and the the 0{a^) have been estimated [1381. ^o the axial anomaly 
of the singlet axial vector current, AS is intrinsically Q^-dependent. Depending on 
the factorization scheme applied [139| this results either in a scale-dependence of the 
sea quark polarization or in an extra contribution to the Ellis- Jaffe sum rule, involving 
Ag = [(7+ {x)—g^ (x)] dx, the gluonic equivalent of the quark distribution moments. 
Both formulations are equivalent. 

Higher twist effects in the dependence of Fi will not be considered here. 

Evaluation of the Fi requires knowledge of gi in the entire interval from to 1. 
Measurements cover a limited kinematic range and thus extrapolations of 171 to and 
1 are necessary. The latter is not critical since gi — s-O at x —^l but the former is 
a considerable problem since gi increases as x decreases and its behaviour at low x 
is theoretically not understood, see sect. 7. 4. Therefore results on Fi depend on the 
assumptions made in the x — s-O extrapolation. Both SMC and SLAC experiments 
assume the Regge like behaviour of gi, i.e. that at x -^0, gi behaves like x", 0< a < 



t Several analyses 1 131 1 incorporating NLO corrections 1 13C ] have appeared after this Workshop. 
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0.5. A value a—0 was chosen and gi was fitted to the two data points at lowest x, 
allowing for variation of this behaviour within the Regge model. This approach might 
however be inconsistent with QCD which predicts a faster rise of gi at low x. 



7.2. Experiments 



New generation polarised electroproduction experiments are listed in table 7.1. 

The experiment of the Spin Muon Collaboration (SMC) at CERN uses a naturally 
polarised muon beam and a cryogenic, solid state target. Experiments E142 - E155 
at SLAC use an electron beam and liquid (solid) cryogenic targets. The HERMES 
experiment at DESY uses an electron beam from the HERA collider and internal 
gas targets. The scattered muon spectrometers in the SMC and SLAC experiments 
have been used (with little change) in DIS experiments preceeding the polarised 
programme, contrary to the HERMES apparatus. Electron and muon measurements 
are complementary: the former offers very high beam intensities and thus statistics 
but its kinematic acceptance is limited to low values of and moderate values of x, 
the latter extends to higher and down to low values of x (an important aspect in 
the study of sum rules) but due to limited muon intensities the data taking time has 
to be long to ensure a satisfactory statistics. 



Table 7.1. New generation experiments on polarised deep inelastic charged lepton — 
nucleon scattering. The las t co lumn shows references to the principal physics results 
obtained until now, (from [140|). 



Experiment 


Beam 


Year 


Beam energy (GeV) 


Target 


References 


SMC 


M+ 


1992-5 
1993 
1996 


100,190 
190 
190 


C4D9OD 
C4H9OH 
NH3 


1 
1 


141, 
142, 


142 
144 


1 
1 


E142 
E143 
E154 
E155 


e" 


1992 
1993 
1995 
1996 


19.4 -25.5 
29.1 

50 
50 


^He 
NH3, ND3 

^He 
NH3, ND3 




[1451 




1 


i4e. 


147 


1 


HERMES 




1995- 


30-35 


H, D, ^He 





The lowest x in the results published by the SMC is about 10~^ and corresponds to 

about 1 GeV^. In the course of analysis are events having lower Q"^ and reaching x 
values of 10~^. A special trigger has been set up recently to extend the measurements 
down to x=10^^, at the expense of lowering to 0.01 GeV^. The upper limit of 

in the SMC is about 100 GeV^. The SLAC experiments' acceptance extends from 
X about 0.01 at ^ 1 GeV^ up ^ iq QeV^ at x ~0.7. As in all fixed-target 
experiments, the low values of x in the SMC and SLAC are correlated with low Q^. 

The cross section asymmetry measured in the polarised lepton - polarised nucleon 
experiments, A^xp, is related to the asymmetries defined in eq.( 
where Pt^Pb denote the target and beam polarisations and /, the target dilution 
factor, accounts for the fact that only a fraction of nucleons is polarised. The beam 
polarisation at the SMC has been measured with a purpose-built polarimeter, using 
two independent methods: polarised jj-e scattering and an analysis of the energy 
spectrum of electrons coming from the muon decay. The result is P^ — — 0.790±0.025 



7. A) hy Ae^p = fPtPb A 
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at 190 GeV beam energy. The target, subdivided into two cells polarised in opposite 
directions, was typically polarised up to 50% for the deuteron and 85% for the 
proton target. The target spin directions were reversed 5 times a day. Polarisation 
of the SLAG electron beam reached 86% in the E143 and was randomly reversed. 
Polarisation of the targets reached 80% for the proton and 25% for the deuteron one 
in E143. Systematic uncertainties in the SMC and SLAG experiments are similar. 



1.3. Results of the measurements and spin structure of the nucleon 



Cross section asymmetries Ai and spin dependent structure functions gi have been 
measured for the proton and deuteron targets by the SMC, |141, 142, 143| , 144 1 and by 
the E143, |146|, |147| . Information on the neutron has been evaluated from the data on 



^He (E142, ||145[) and from the data on the proton and deuteron (SMC, ||l4l|, |l42(|). 
All data sets are in a very good mutual agreement even if Ai, extracted from data 
covering different intervals, has been assumed to be independent. 

Results on gi for proton, deuteron and neutron are shown in Figs. 7.1, 7.2. Here 
gi = 2gf/(l — 1.5ujd) — g\ where lud ^ 0.05 is the probability of the D-state of the 
deuteron. Conversion of Ai to gi, which was made under an assumption that Ai scales, 
needs only information on the structure function Fi or, equivalently, F2 and R (cf. 
eq.( |7.6| )). The NMG parametrisation of ^2(2;, Q^) [Q and the SLAC parametrisation 
of R(x, Q2) pg[] have been used by both SMC and SLAC. 



2.5 



gf(x) 



1.5 





-0.5 



T E80/E130 

▲ EMC 

• SMC 

□ £143 

Q^= 5 GeV^ 



Figure 7.1. The spin dependent structure function gi{x) of the proton at Q^=5 
GeV'^. the EMC data were reevaluated using the same F2 and R parametrisations 
as for the SMC and E143 data. Error bars are s tatisticah the shaded area marks the 
SMC systematic errors. Figure taken from |l40| . 

The behaviour of the is different from that of gf and gi, especially at low x. 
This should be contrasted with the unpolarised case where a small difference between 
proton and neutron structure functions can be explained by nuclear shadowing in 



the deute ron. Measurements of the asymmetry, A2, for the proton |144 and the 



deuteron [148| showed that this function is significantly smaller than the bound ^fR 
and consistent with zero. 

The Bjorken sum rule seems to be fulfilled by the above data at the 10% level: its 
value measured by the SMC at =10 GeV^ is Tl-V'l =0.199±0.038, to be compared 
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0.5 
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a 4 ^» 



b) 



= 5 GeV" 

• SMC 

□ E143 
A E142 



10 



Figure 7.2. The spin dependent structure functions (a) gf(x) and (b) g"{x), as a 
function of x at Qg=5 GeV^. error bars a re st atistical: the shaded area marks the 
SMC systematic errors. Figure taken from ll42|. 



with the QCD prediction (four flavours, corrections up to ag): 0.187±0.003. The EUis- 
Jaffe sum rule is not confirmed by the data, the pattern of disagreement being similar 
in the proton and deuteron results. At = 5 GeV^, the combined data of the 
EMC, SMC and SLAG give for rP,Tf and respectively: 0.125±0.009, 0.041±0.005 
and -0.037±0.008 as compared to the predicted 0.167±0.005, 0.070±0.004 and 
— 0.015±0.005. A most straightforward explanation of this violation may be a non- 
zero polarisation of the strange sea. Results of the new SMC proton data analysis, 
with extended kinematic coverage and a NLO QCD analysis to evolve the measured 
gi{x, Q^) to a common value of Q^, confirm all the above conclusions. 
The nucleon spin, Sz — ^, can be decomposed as follows 

Sz = ^AY. + Ag + Lz (7.10) 

where Lz is angular momentum due to the partons. All the data sets (except perhaps 
E142), evaluated with a consistent treatment of the QCD corrections at a common 
GeV^ and under the assumption that the flavour SU(3) is exact, show that the 
total quark contribution to the nucleon spin is small, AS ~0.2, and that the strange 
sea is indeed polarised: As ^ —0.1. SU(3)-breaking can decrease As but leaves AS 
unchanged. Choosing a factorization scheme in which the quarks polarisation is scale 
independent, a dependent gluonic contribution appears in the Ellis- Jaffe sum rule 



as a result of the anomalous dimension of the singlet axial vector current |139|. Then 
the Ellis- Jaffe assumption of As=0 implies that at GeV^, A^ 3 is needed to 

restore the sum rule. 

Finally we note the first measurements of the semi-inclusive spin asymmetries for 
positively and negatively charged hadrons in the polarised muon-proton and muon- 



deuteron scattering in the SMC experiment |149|. The x dependence of the spin 
distributions for the up and down valence quarks and for the non-strange sea quarks 
has been determined. The moments of the quark spin distributions were obtained to 
be: Au„=1.01±0.24, Ady = — 0.57±0.25; moments for the non-strange sea quarks are 
consistent with zero over the whole measured range of x. 
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7.^. Theory predictions for gi at small x 

Observation of a difference between the proton and deuteron spin structure functions 
g\''^{x, Q^) at small x, Figs. 7.1, 7.2, indicates a sizeable non-singlet contribution to gi 
in this region. The dominant behaviour of g" '*' can be estimated from a resummation 
of the leading QCD corrections. The conventional Altarelli-Parisi approach only 
resums terms of the form (as In ln(l/a;))'^ and yields 



5r(-,Q^)-exp(^^^ln^lnij. (7.11) 
Terms proportional to (as ln^(l/a;))'^, which arise due to the presence of two fermions 



in the t-channel |150| are, although more singular at small cc, not included. 



Ryskin outlined a calculation of the small x behaviour of , which resums these 



contributions |151]. This calculation not only includes the enlarged (compared to the 
Altarelli-Parisi approach) phase space of ladder diagrams which leads to the double 
logarithmic contributions, but also substantially differs from similar resummations 
performed for unpolarised structure functions in that both ladder and non-ladder 
diagrams contribute to the most singular part of the amplitude. Using the fact that 
the non-ladder gluon has to be softer than all gluons within the ladder, these non- 
ladder contributions can be resummed. This procedure results in a nonlinear infrared 



evolution equation |152| for the quark amplitudes of odd signature, which mixes color- 
singlet and color-octet quark amplitudes. By means of a Mellin transformation, this 
equation can be solved, yielding 

gu.s. _ ^-4-) ^i^i^ _ iM^/2asCF/TT ^ 0.4 (for as = 0.18). 

(7.12) 

The non-ladder (bremsstrahlung) contributions give a relatively small contribution 
(4%) to wiri; if only ladder diagrams are retained then w^Ti = uj^}. = ^^asCp j-n for 
fixed coupling. When running coupling effects are included in the double logarithmic 
resummation generated by the ladder diagrams then the effective slope is reduced to 
0.2 - 0.3 [HI. 

First preliminary results of a similar calculation for the singlet contribution to 



g\ [x] at small x were presented [ |154[ , indicating that 

with ^^7^ ~ 3 ^^72 > 1 (for as > 0.12). (7.13) 

As the above result contradicts conventional Regge-type extrapolations of g\^^^ in the 
EUis-Jaffe sum rule, it shows the need for extending the measurements of (7i(a;,Q^) 
down to the lowest possible values of x. 



7. 5. Prospects for the future 

Understanding of the polarised structure functions has improved dramatically in the 
recent years, thanks to the EMC, SMC and SLAC measuremets. Several questions 
however remained unanswered. Among them is the low x behaviour of gi, its 
evolution, the gluon polarisation and flavour decomposition of polarised parton 
distribution. The HERMES experiment, recently starting at HERA, using a polarised 
electron beam and a polarised internal gas target will especially address the last 
question from a presently unique reconstruction of the hadronic final state. To answer 
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the remaining questions a new generation of experiments, e.g. at the HERA collider, 
is needed. Prospects of spin physics at HERA were discussed at a workshop at 
DESY-Zeuthen in August 1995. A polarised deep inelastic programme at HERA 
could allow measurements over an extended kinematic range, including low x and 
high . Polarisation of the proton beam is technically much more complicated than 
polarisation of the electron beam, as the proton beam does not polarise naturally 
[155|. Construction of the polarised proton beams of energy up to 250 GeV in the 
RHIC collider rings has already been approved, a helpful step for HERA. Various 
suggestions for dedicated measurements of Ag^XjQ"^), including the HMC /CHEOPS 
project at CERN, were also discussed at the Zeuthen workshop [156|. 



8. Conclusions 



We believe that this report gives a fair summary of our present understanding of the 
structure of the proton, as it is seen in deep-inelastic lepton-proton scattering. Much 
— but not all — of the recent experimental information comes from HERA. We have 
summarized the HERA and fixed-target data in the first part of this Report. 

The most puzzling phenomenon continues to be the rise of F2 at small x, which 
seems to persist even at low values Q^, and its connection with the soft Pomeron 
seen in photoproduction at = 0. So far there is no real problem in describing the 
inelastic data within the CLAP evolution scheme, i.e. it is possible to find starting 
distributions for quarks and gluons which evolve according to the (next-to-leading) 
QCD evolution equations and successfully describe the data. This alone, however, 
does not explain the observed rise: either one has to ascribe the origin of the rise 
to the input (gluon) distribution (the BFKL scenario), or one has to allow a rather 
large range of the CLAP evolution. In the latter case, the GLAP evolution by 
itself leads to a rise a small x. But in order to accommodate the observed rise at 
low the QCD evolution has to start at an even lower scale where - on rather 
general grounds - one would have been hesitant to use a leading-twist perturbative 
QCD framework. A further ambiguity is that, as yet, the "BFKL scenario" is based 
only on the resummation of leading log(l/a;) terms. A quantitative description needs 
the computation of sub- leading effects. 

The discussions of the Working Group, as described in this summary, reflect the 
agreement among the participants that we need a thorough understanding of the 
validity of GLAP and BFKL descriptions (and a description which unifies both) in 
the low- a; and low-Q^ region, as well as the interpolation between deep inelastic 
scattering and photoproduction, before a deeper understanding of the interface 
between perturbative and non-perturbative physics can be reached. It seems clear 
that to gain further insight we must study more than just the fully inclusive total 
cross section as measured by F2, i.e. we need to look more closely into the final state 
of deep-inelastic scattering. 

The discussion of the spin structure of the proton also seems to be developing in 
the small-x direction, somewhat analogous to the unpolarized case. Namely, the low- a; 
data points for both the proton and the neutron spin structure functions 171 indicate 
the possibility of a rise at small x. At the same time QCD calculations, based on 
leading and next-to-leading order Altarelli-Parisi splitting functions or, more recently, 
using the double-logarithmic approximation, predict at strong rise of gi at small x. 
So there is little doubt that small- a; physics is becoming a field of particular interest 
also in the proton spin community. 
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